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ABSTRACT 


Comparisons  of  theory  and  experiment  have  shown  that  the  vector 
electromagnetic  theory  of  scattering  from  microrough  surfaces  is  largely  in  agreement 
with  measurements.  Data  were  taken  at  several  wavelengths,  for  several  angles  of 
incidence,  and  for  several  combinations  of  polarization.  The  agreement  for  the  most 
part  is  within  experimental  error.  Surface  data  were  obtained  on  the  samples  by 
stylus  profilometry  and  interferometric  probes.  The  data  were  used  to  conceive  and 
design  an  optical  shop  scatterometer.  using  a  helium-neon  laser  and  ten  silicon 
detectors.  The  instrument  was  shown  to  be  capable  of  measuring  mirror  samples  of 
reasonable  roughness,  and  scaling  to  other  wavelengths  could  then  be  accomplished. 
Thus,  a  simple  instrument  could  be  used  in  the  optical  shop  to  characterize  mirrors 
up  to  2  m  in  diameter.  Very  good  mirrors  could  not  be  scaled.  The  rollover  point 
in  the  curves  was  too  close  to  the  specular  direction;  values  of  the  surface  roughness 
and  correlation  length  could  not  be  determined  independently. 

A  second  instrument,  operating  at  about  10  /im.  was  designed  and  a  laboratory 
model  was  constructed.  A  noise  equivalent  BRDF  of  lO^sr*1  at  2°  from  specular 
was  demonstrated  with  an  instrument  that  has  1-m-long  arms.  Thus,  it  was  shown 
that  a  device  consisting  of  a  C02  laser  and  pyroelectric  detectors  can  be  constructed 
that  is  capable  of  characterizing  mirrors  of  very  good  quality. 

Appendices  give  design  details  and  instructions  for  use. 
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INTRODUCTION 


This  is  the  final  report  on  a  study  by  the  University  of  Arizona  for  the  Rome 
Air  Develoment  Center  at  Griffiss  Air  Force  Base.  The  study  had  as  its  primary 
goal  the  development  of  a  scatterometer  that  could  be  used  in  an  optical  shop  on 
finished  mirrors  with  diameters  as  large  as  two  meters.  One  of  the  measurement 
goals  was  a  bidirectional  reflectance  distribution  function  (BRDF)  of  10-3sr“*  one 
degree  from  specular  at  10.6  fim.  The  details  of  the  work  statement  are  given  in 
Appendix  A. 

Theoretical  developments  were  to  be  carried  out  to  do  the  appropriate 
wavelength  and  angle  scaling.  They  were  to  be  accompanied  by  appropriate 
measurements. 


PRINCIPAL  RESULTS 

Two  versions  of  a  scatterometer  were  constructed  using  a  helium-neon  laser, 
chopper,  silicon  detectors,  and  appropriate  electronics.  One  version  simply  used  ten 
detectors,  properly  spaced  to  obtain  the  scatter  measurements.  This  instrument  is 
described  in  a  published  paper  that  is  reproduced  as  Appendix  B.  The  instrument 
was  demonstrated  at  RADC  in  1986.  Instructions  for  its  use  are  given  in 
Appendix  C. 

One  shortcoming  of  that  instrument  was  its  limitation  in  obtaining  angularly 
resolved  data  near  the  specular  direction.  Therefore,  a  second  instrument  was 
delivered  and  demonstrated  in  1987.  It  used  an  array  of  silicon  detectors  in  place  of 
the  single  detector  closest  to  the  specular  reflection  direction.  This  instrument  was 
described  in  another  paper.  A  copy  is  attached  as  Appendix  D.  There  was  still  a 
problem  resulting  from  instrument  background. 

A  third  instrument  was  therefore  designed.  It  was  based  on  a  C02  laser  and 
ten  pyroelectric  detectors.  The  same  geometry,  a  chopper,  and  appropriate  electronics 
were  designed.  Because  it  was  imperative  that  good  performance  near  specular  be 
achieved,  we  carried  out  laboratory  experiments.  We  obtained  a  noise  equivalent 
BRDF  of  10-4sr-i  at  2°  from  specular  at  10.6  /um.  A  presently  unpublished  article 
describing  this  is  enclosed  as  Appendix  E. 
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TECHNICAL  DISCUSSION 


This  section  provides  some  of  the  results  of  the  theory  of  scattering  so  that  the 
instrument  designs  and  limitations  may  be  better  understood. 

Some  Scaling  Laws 

Several  authors  have  published  theoretical  results  of  scattering  from  surfaces  with 
roughness  less  than  the  wavelength  of  the  radiation.1-15  The  basic  result  is 

Pb  ~  BRDF 
Pb  -  F(X)F(0)F(s) 

F(X)  -  wavelength  factor 
F(0)  -  angle  factor 
F(s)  -  surface  factor 

BRDF  -  scattered  radiance  divided  by  incident  irradiance. 

The  BRDF  may  be  written  as  the  product  of  three  factors:  wavelength,  geometry, 
and  surface.  The  first  of  these  is  either  1A3  or  Ifk*  depending  on  the  details  of  the 
model.  The  second  is  a  function  of  the  dielectric  constant  and  angles  of  incidence 
and  scatter.  It  is  a  relatively  weak  function,  and  it  may  be  ignored  in  the  ensuing 
analysis.  The  third  factor  relates  directly  to  the  surface.  It  is.  in  fact,  the  variance 
spectrum  of  the  surface-height  distribution.  For  a  random  surface  with  exponential 
autocorrelation,  it  can  be  shown  to  have  one  of  two  spectra: 


where 


Pb 


k3F  (6)or2i  k4F(9)g2f2 

l+P2*2  (1+p2!2)3/2 


k  -  2nfk 

g  -  rms  roughness 
I  -  correlation  length 
P  -  k(0-0o) 

-  k(sin<?-sin0o) 

9  m  scatter  angle 
90  -  specular  angle. 


The  first  relates  to  a  one-dimensional  transform;  the  second  to  a  two-dimensional 
version.  Most  data  fit  the  first  form. 


We  are  interested  in  how  these  scale  with  wavelength  and  angle.  We  ignore 
F (0),  that  is,  we  let  F(d)  be  1,  and  investigate  the  linear  asymptotes  of  these 
expressions,  for  p2f2«l  and  p2f2))  1 .  In  the  first  case 

Ph  -  kVi  or  ^  -  k4cr2i2  . 

The  BRDF  is  a  straight,  horizontal  line.  The  two-dimensional  version  is  2irifk  times 
the  one-dimensional  one.  We  shall  see  that  this  is  often  not  significant.  In  the  other 
region 


BRDF 

BRDF 


2k  3<r2i  2  a2 

ffp2i2  "  ff\i(/3-/30)2 

2k*g2i2  2j2 

p3!3  “  XiOS-/J0)3  • 


Finally,  the  asymptotes  for  both  forms  intersect  at  p  -  1/1.  Several  things  arise  from 
these  theoretical  results:  the  asymptotic  curves  are  straight,  horizontal  lines  at  small 
angles  and  they  are  slopes  of  either  -2  or  -3  at  larger  angles  when  plotted 
logarithmically.  The  intersection  is  at  p  -  1/1. 

Validity  of  Scaling  Theories 

Several  theories  were  investigated,  and  appropriate  data  were  taken  to  test  these 
theories.  This  information  has  been  published  in  the  Journal  of  the  Optical  Society 
of  America,  and  a  reprint  is  included  here  as  Appendix  F.  The  results  are.  in 
essence,  for  specular  surfaces,  that  the  vector  electromagnetic  theory  is  valid;  the 
scalar  theory  is  not.  For  measurements  in  the  plane  of  incidence,  the  one¬ 
dimensional  transform  applies;  the  two-dimensional  version  does  not. 

Appendix  G  is  a  reprint  from  an  SPIE  Proceedings  that  also  discusses  the 
inference  of  surface  parameters  by  various  techniques.  The  one  we  prefer  is  fitting 
the  shape  of  the  curve  with  both  the  correlation  length  and  the  rms  roughness. 

Measurements  of  Surface  Parameters  from  BRDF  Values 

Two  values  can  be  inferred  from  the  BRDF  vs  0  -  0O  curves,  like  these,  the  rms 
roughness  g  and  the  autocorrelation  length  1.  Thus  we  need  at  least  two 
measurements.  If  we  do  it  on  the  slope,  we  have 


The  subscripts  1  and  2  represent  measurement  values.  We  assume,  as  is  the  custom, 
that  these  were  both  measured  at  the  same  wavelength  and  that  neither  a  nor  f 
changed.  Then  we  have  the  simple  relationship 


*>,  0S-/3o)l 

Pb,  '  owT  ‘ 

They  are  linearly  related.  This  means  that  to  infer  I,  and  therefore  or,  we  must  at 
least  approach  the  region  where  p  -  1/1.  A  typical  correlation  length  for  a  good 
surface  is  10  jtim.  Then  at  0.6328  /xm . 

/J-/30  -  \/2iri  -  0.01  . 

If  we  also  use  an  rms  roughness  of  10  nm  (0.01  /xm).  we  can  find  the  BRDF  at  this 
spatial  frequency: 


BRDF  -  k3or2IA 

-  [oik]!(0-01>’10/2 

as  0.49  sr*‘  . 

This  is  a  difficult  measurement  to  achieve.  The  specification  of  10-3sr-1  at  1° 
means  that  the  product  of  the  rms  roughness  and  the  correlation  length  must  be  about 
10"6  /xm3.  This  is  equivalent  to  a  mirror  with  a  1-nm  roughness  and  a  l-/xm 
correlation  length,  or  a  0.3-nm  roughness  and  10-/xm  correlation  length,  very  good 
mirrors  indeed.  In  fact,  on  some  of  the  better  mirrors  the  intersection  could  not  be 
determined. 

The  intersection  point  is  wavelength  dependent.  In  fact,  it  is  linearly  dependent 
upon  the  wavelength.  Thus,  in  moving  from  0.63  /xm  to  10.6  /xm,  the  0  -  /3„  value  of 
the  intersection  should  move  from  0.01  to  about  0.17.  However,  the  flatband  BRDF 
will  be  reduced  by  a  factor  of  almost  five  thousand.  The  measurement  is  better 
made  at  10  /xm,  which  is  why  we  developed  the  infrared  instrument.  The  infrared 
instrument,  as  is  shown  in  Appendix  E,  would  measure  a  BRDF  of  10_3sr_1  with 
about  10%  error  as  close  to  specular  as  2°  Q3  -  jS0  -  0.035).  Lower  BRDF  values 
would  be  possible  at  larger  angles.  The  noise  equivalent  BRDF  of  about  6x1 0-6  sr"1 
would  be  obtained  at  about  6°  (/J  -  ft,  -  0.1).  This  is  well  below  the  indicated 
rolloff  values.  For  instance,  a  sample  with  a  correlation  length  of  10  /im  and  an  rms 
roughness  of  1  nm  has  roll-off  values  of  /3  -  0O  -  0.16  and  a  corresponding  BRDF 


value  of  1.6x10"*  sr"1.  Unfortunately,  there  were  insufficient  funds  for  the 
construction  of  a  final  infrared  unit;  we  were  only  able  to  show  feasibility. 

Non-Random  Surfaces 

The  vector  version  of  the  electromagnetic  theory  in  its  basic  form,  which  we  have 
demonstrated  to  be  the  most  accurate  representation  of  smooth  surfaces,  requires  the 
spectrum  of  the  surface-height  distribution.  Diamond-turned  samples  do  not  have 
randomly  polished  surfaces.  They  may  be  characterized  by  the  very  fine  grooves  of 
the  diamond  tool  as  it  advances.  They  are,  in  fact,  weak  gratings  with  a  spatial 
frequency  of  p  -  k(/S  -  /S0).  Thus  (radian)  spatial  frequency  is  related  to  the  spatial 
period  s. 


p  -  k(jS  -  £0)  -  2?r/s  . 

Then 

s  -  X/0S-/3o) 

/S-/S0  -  X/s 

for  the  equivalent  correlation  length.  Thus,  for  example,  for  a  wavelength  of 
0.63  nm  and  screw  advance  of  1  fim,  there  will  be  a  peak  at  p  -  /3„  -  0.6.  This  is 
well  into  the  region  of  effective  measurement.  For  the  peak  to  be  at  less  than  0.1, 
the  advance  would  have  to  be  more  than  6.3  fim.  Some  of  our  BRDF  data  have 
been  used  to  calculate  this  screw  advance  with  good  agreement. 

CONCLUSIONS 

The  vector  electromagnetic  theory  is  a  good  model  for  microrough  surfaces.  It 
incorporates  three  factors,  the  most  important  of  which  is  the  surface-height  variance 
spectrum.  It  applies  to  both  random  and  to  diamond-turned  surfaces  in  the  plane  of 
incidence.  An  instrument  has  been  designed  and  constructed  that  can  measure  the 
BRDF  of  large  samples  at  various  points  at  a  wavelength  of  0.6328  fim.  It  can 
measure  a  BRDF  of  approximately  10"s  sr-1  as  close  as  about  0  -  0O  -  0.01.  Thus,  it 
should  be  good  for  samples  that  have  a  correlation  length  of  about  10  fim  and  an  rms 
roughness  of  about  0.4  nm.  Samples  better  than  this  can  be  measured  but  not 
scaled.  An  infrared  scatterometer  is  feasible.  A  laboratory  demonstration,  using  the 
geometry  and  components  that  would  be  incorporated  in  a  final  model,  showed  that  a 
noise-equivalent  BRDF  of  10"*  sr"1  was  obtained  at  an  angle  of  2°  from  the  specular. 


BIBLIOGRAPHY  OF  VECTOR  SCATTERING  THEORY 


1.  Elson,  J.  M.,  "Light  scattering  from  semi-infinite  media  for  non-normal 
incidence,"  Phys.  Rev.  B  12,  2541  (1975). 

2.  Church,  E.  L.,  H.  A.  Jenkinson,  and  J.  M.  Zavada,  "Measurement  of  the  finish 
of  diamond-turned  metal  surfaces  by  differential  light  scattering,"  Opt.  Eng.  16,  360 
(1977). 

3.  Elson,  J.  M.  and  R.  H.  Ritchie,  "Diffuse  scattering  and  surface-plasmon 
generation  by  photons  at  a  rough  dielectric  surface,"  Physicus  Status  Solidi  62,  461 
(1974). 

4.  Maradudin,  A.  A.  and  W.  Zierau,  "Effects  of  surface  roughness  on  the  surface- 
polariton  dispersion  relation,"  Phys.  Rev.  B  14,  484  (1976). 

5.  Maradudin.  A.  A.  and  D.  L.  Mills,  "Scattering  and  absorption  of  electromagnetic 
radiation  by  a  semi-infinite  medium  in  the  presence  of  surface  roughness,"  Phys. 
Rev.  B  11.  11392  (1975). 

6.  Kroger,  E.  and  E.  Kretschmann,  "Scattering  by  slightly  rough  surfaces  or  thin 
films  including  plasma  resonance  emission,"  Zeitschrift  fur  Physik,  237,  1  (1970). 

7.  Hill,  N.  R„  "Integral-equation  perturbative  approach  to  optical  scattering  from 
rough  surfaces,"  Phys.  Rev.  B  24,  7112  (1981). 

8.  Marvin,  A..  F.  Toigo,  and  V.  Celli,  "Light  scattering  from  rough  surfaces: 

General  incidence  angle  and  polarization."  Phys.  Rev.  B  11,  2777  (1975). 

9.  Church,  E.  L.,  H.  A.  Jenkinson,  and  J.  M.  Zavada,  "Relationship  between  surface 
scattering  and  microtopographic  features,"  Opt.  Eng.  18,  125  (1979). 

10.  Church,  E.  L.  and  J.  M.  Zavad,  "Residual  roughness  of  diamond  turned  optics," 
Appl.  Opt.  14,  1788  (1975). 

11.  Rice,  S.  O.,  "Reflection  of  electromagnetic  waves  from  slightly  rough  surfaces," 
Communications  in  Pure  and  Applied  Mathematics  4,  351  (1951). 

12.  Bennett,  M.  E.  and  J.  O.  Porteus,  "Relation  between  surface  roughness  and 
specular  reflectance  at  normal  incidence,"  J.  Opt.  Soc.  Am.  51,  123  (1961). 

13.  Elson,  J.  M.  and  H.  E.  Bennett,  "Relation  between  the  angular  dependence  of 
scattering  and  the  statistical  properties  of  optical  surfaces,"  J.  Opt.  Soc.  Am.  69.  31 
(1979). 

14.  Crowell,  J.  and  R.  H.  Ritchie,  "Surface-plasmon  effect  in  the  reflectance  of  a 
metal,"  J.  Opt.  Soc.  Am.  60,  794  (1970). 

15.  Toigo,  F.,  A.  Marvin,  V.  Celli,  and  N.  R.  Hill,  "Optical  properties  of  rough 
surfaces:  General  theory  and  the  small  roughness  limit,"  Phys.  Rev.  B  15,  5618 
(1977). 


APPENDIX  A 


STATEMENT  OF  WORK 

The  University  of  Arizona,  in  the  development  of  the  bidirectional  reflectance 
distribution  function  (BRDF)  scatterometer,  shall  accomplish  the  following  tasks: 

a.  Compare  established  scatter  theories  to  be  able  to  predict  a  BRDF  of 
10“3  sr-1  at  1°  from  specular  at  10.6  fim.  The  modeling  of  these  theories 
shall  be  such  that  both  the  wavelength  scaling  to  10.6  nm  and  angular 
scaling  to  1°  are  achievable. 

b.  Adapt  the  Arizona  Scatterometer  (AZSCAT)  to  measure  steeply  sloped  and 
high-order  aspheric  ssurfaces  at  a  minimum  BRDF  of  lO^sr"1  at  1°  from 
specular  at  10.6  fim. 

c.  Measure  appropriate  specimens  to  adequately  verify  modeling  of  the  scatter 
theory  and  to  assess  the  adaptation  of  AZSCAT  to  the  parameters  specified 
in  paragraph  b  above. 

d.  Design,  fabricate,  test  and  deliver  a  protable  scatterometer  that  is  capable  of 
measuring  or  using  appropriately  scaled  data  to  infer  BRDF  of  10_3sr_1  or 
less  at  1°  from  specular  at  10.6  fim.  The  instrument  shall  be  capable  of 
measuring  BRDF  on  a  surface  with  F/l  or  less  and  a  high-order  aspheric. 

e.  Verify  the  utility  of  the  instrument  to  the  design  parameters  using  samples 
supplied  by  the  Contractor  and,  if  desired  by  the  Government,  those  samples 
directed  by  the  Contracting  Officer. 

f.  Provide  familiarization  to  a  maximum  of  ten  Government  selected  personnel 
in  the  application  of  data  and  the  instrument  to  analyze  the  stray  light 
properties  of  developmental  or  operational  systems,  at  no  more  than  three 
instructional  sessions. 
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APPENDIX  B 


A  Portable  Scatterometer  for  Optical  Shop  ties 

Hill lam  £,.  Wolfe,  Kathleen  Magee  and  Douglas  H.  Wolfe 

Optical  Science#  Center,  University  of  Arizona 
Tucson,  Arizona 


In t roduc t Ion 

The  scattering  of  optical  components  has  become  a  more  and  aore  critical  Issue 
In  modern  sensor  systems.  This  Is  a  result  of  Improved  performance  capabilities  on 
the  part  of  detectors  and  of  optical  elements  themselves,  as  well  as  an  Increased 
awareness  of  the  effects  of  scatter  on  system  performance.  There  are  several 
programs  that  evaluate  scatter  in  complicated  optical  trains,  the  best  known  of 
which  are  APART  and  GUERAP  (Arizona  Paraxial  Analysis  of  Radiative  Transfer  and 
General  Unwanted  Energy  Radiation  Analysis  Program).  These  programs  require 
component  scatter  data  as  Inputs  f o  .•  the  analysis  and  even  the  beginning  design 
approaches.  These  data  have  generally  been  generated  with  laboratory  Instruments 
that  were  essentially  gonlometrlc  radiometers.  The  data  have  been  difficult  to 
obtain  and  usually  obtained  for  small  representative  samples  of  the  real  surfaces. 
If,  In  the  process  of  manufacture  of  the  parts  a  measurement  of  the  scatter  was 
made.  It  was  done  with  a  relatively  crude  Jerry-rigged  apparatus.  Sometimes  the 
measurements  made  this  way  were  quite  accurate,  but  they  were  never  easy. 

This  paper  describes  the  design,  construction  and  performance  of  a  portable 
apparatus  that  can  make  scattering  measurements  on  samples  as  large  as  2  meters  In 
diameter  at  0.6328um  for  a  45-degree  angle  of  incidence  and  a  reasonably  large 
range  of  scatter  angles. 

Theoretical  background 

The  scattering  of  mirror  surfaces  is  usually  described  In  terms  of  the 
Bidirectional  Reflectance  Distribution  Function  (BRDF).  Transmissive  samples  may  be 
characterized  by  their  values  of  Bidirectional  Transmission  Distribution  Function 
( BTDF )  and  the  general  function  can  be  the  BSDF .  where  the  S  is  for  scattering.  The 
Instrumentation  to  be  described  here  Is  a  BRDF  measurement  device,  meant  for  mirror 
samples,  although  it  can  be  used  equally  well  for  diffuse  reflectors.  In  Its 
present  configuration  it  cannot  be  used  for  the  measurement  of  transparent-sample 
BTDF  values . 

Most  of  the  data  obtained  on  the  BRDF  values  of  both  specular  and  diffuse 
samples  have  beer,  presented  In  a  log-log  format  in  which  the  BRDF  is  plotted  as  a 
function  of  the  surface  spatial  frequency  for  a  fixed  wavelength.  The  coordinate 
is  commonly  referred  to  as  "beta-beta  naught"  or  "delta  beta."  An  example  is  shown 
in  Fig.  1;  it  is  the  BRDF  of  a  mirror  measured  at  10.6um.  The  values  range  from 
almost  one  to  about  one-millionth  in  units  or  reciprocal  steradians.  This  of 
course  is  the  reason  for  the  logarithmic  ordinate.  The  abscissa  ranges  from  0.01  to 
2.  The  maximum  possible  value  for  the  abscissa  is  2,  while  the  minimum  theoretical 
value  is  zero.  The  minimum  practical  value  for  such  measurements  is  about  0.018  or 
about  one  degree.  Again,  the  range  is  more  than  an  order  of  magnitude,  and  the 
small  values  are  Just  as  important  (or  more  important)  than  the  large  ones. 

A  diffuse  sample  will,  of  course,  require  the  same  range  for  the  independent 
value,  but  have  a  much  smaller  range  of  BRDF  values.  An  example  of  a  good,  diffuse 
black  is  shown  in  Figure  2.  Nothing  is  lost  by  presenting  this  information  on  a 
logarithmic  scale. 

Vector  electromagnetic  theory  has  been  used  to  model  the  scattering  from 
surfaces  like  these,  which  have  roughness  values  much  less  than  the  wavelength  of 
the  Incident  light.  The  theory  predicts  that  the  BRDF  can  be  written  as  the 
product  of  three  factors,  a  wavelength  factor,  a  geometry  factor  and  a  suface 
factor: 

GBRDF-KGW 

For  measurements  in  the  plane  of  Incidence  the  wavelength  factor  is  K  cubed, 
where  k  Is  2  pi  divided  by  the  wavelength.  The  geometry  factor  is  a  function  of 
olarizatior.:  it  is  a  smooth  function  of  delta  beta  at  least  until  the  extremes  of 
he  measurement  for  ss  and  pp  polarizations,  and  it  is  zero  for  the  crossed 
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polarizations  in  the  plane  at  incidence.  These  factors  are  plotted  in  Fig.  3.  The 
unpolarized  version  is  the  average  of  these  two.  The  point  is  that  these  are 
smooth  functions.  The  surface  factor  W  is  the  power  spectrum  of  the  surface  heigh; 
distribution.  For  most  mirrors  made  in  the  classical  way  this  is  a  simple 
Lorentian  that  can  be  characterized  by  two  asymptotes,  one  that  is  a  horizontal 
straight  line,  and  one  that  is  a  straight  line  with  a  slope  of  minus  two.  Surfaces 
that  have  been  formed  by  diamond  turning  have  superimposed  on  these  curves  a  set  of 
peaks  that  correspond  to  the  advance  of  the  tool. 

The  values  for  scattering  out  of  the  plane  of  incidence  have  not  yet  been 
mode  1 ed . 

Instrumental  ..overview 

Consideration  of  the  data  presented  in  these  figures  indicate  that  the  BRDF  data 
can  be  obtained  by  illuminating  the  samples  with  laser  light  of  the  proper 
wavelength  and  measuring  with  an  array  of  detectors;  it  should  not  be  necessary  to 
have  a  continuous  distribution  of  them.  Data  for  other  angles  of  Incidence  and 
other  wavelengths  can  be  obtained  by  the  scaling  laws  generated  and  developed  from 
the  theory.  These  are  the  basic  concepts  of  the  Instrument  design.  The  questions 
which  arise  are  the  placement  of  the  detectors  and  the  source  with  respect  to  the 
sample,  the  signal  to  noise  ratios,  the  size  of  the  illuminating  spot,  calibration, 
orientation.  data  taking  and  presentation  of  the  results.  We  have  not  tried  to 
design  a  research  instrument,  but  rather  one  that  can  determine  the  properties  of 
most  mirrors  and  diffusers  at  most  wavelengths  for  a  reasonable  range  of  scatter 
ang les . 

The  instrument  consists  of  a  helium-neon  laser,  a  chopper,  a  folding  flat 
mirror,  four  silicon  detectors  with  apertures,  six  with  lenses,  and  amplifiers. 
The  outputs  of  the  amplifiers  are  sampled  with  a  microcomputer  and  displayed  on  its 
screen.  Options  for  printing  the  output  are  therby  provided. 

Geome  t  ry 

The  data  should  be  taken  at  equal  logarithmic  increments  from  a  value  of  0.01  to 
1.  The  values  out  to  2  are  not  really  needed,  and  values  below  0.01  are  Impossible 
to  obtain.  Values  near  0.01  are  indeed  difficult.  It  appears  that  ten  data  points 
across  this  range  are  adequate  because  of  the  smoothness  of  most  of  the  samples. 
Table  1  shows  the  detector  placement  for  these,  including  a  detector  placed  at  the 
specular  position. 

But  what  angle  of  incidence  should  be  used?  If  the  sample  is  illuminated 
normally,  then  the  specular  direction  cannot  be  measured.  If  the  angle  of 
incidence  is  too  large,  then  the  spot  on  the  sample  is  too  large,  although  this  is 
not  a  serious  constraint;  the  laser  spot  is  collimated  into  a  14mm  beam.  The 
incidence  angle  is  really  dictated  by  the  practical  issues  of  how  large  the 
scatteroraeter  can  be  and  how  far  it  must  be  from  the  test  piece.  These 
considerations  must  be  traded  off  with  the  available  range  of  delta  beta  for 
different  angles  of  incidence.  Either  forward  scatter  or  backscatter  or  some  of 
each  may  be  used.  If  the  minimum  practical  incidence  angle  of  20  degrees  is  used 
then  the  range  is  0  to  0.66.  This  is  not  as  large  as  we  would  like.  The  range  is 
smaller  for  70  degrees,  and  the  range  for  10  degrees  is  still  only  0  to  0.83. 
Backscatter lng  seems  more  appropriate.  The  range  then  is  from  0  to  -(1+sinA) .  For 
20  degrees  this  Is  0  to  -1.34;  for  10,-1.17;  for  30, -1.5;  for  45,-1.71.  Any  angle 
of  incidence  gives  a  wide  range  of  delta  beta. 

_P  ra ctl cal.  geometries . 

The  system  should  be  relatively  easy  to  use  in  the  optical  shop  for  mirrors  up 
to  about  2  meters  in  diameter.  We  have  taken  about  60  centimeters  as  the  required 
height  clearance  for  the  system.  we  have  chosen  to  mount  all  the  components  on  an 
aluminium  optical  rail  in  the  shape  of  a  tee.  This  structure  is  supported  by  three 
commercial  tripods.  The  leg  of  the  tee  is  almost  2  meters  long,  while  the  crossbar 
is  about  1.2  meters.  This  leaves,  allowing  for  some  clearances,  about  130cm  for 
the  distance  between  the  input  beam  and  the  specular  return.  For  the  clearance  we 
have  chosen  and  the  length  we  used,  the  angle  of  Incidence  must  be  Just  about  45 
degrees.  We  set  it  at  that  value,  and  then  we  readjusted  the  assumptions  a  little. 

Signal  to  noise  cons lderat ions 

The  noise  equivalent  8R0F  is  given  by  the  following  equation: 
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The  geometry  chosen  for  the  system  results  In  a  distance  from  the  detector  to 
the  spot  of  about  85cm  for  the  specular  direction,  and  this  means  that  to  obtain  a 
solid  angle  0  of  IE-4  the  collecting  area  must  be  0.85  cm.  This  would  be  to  obtain 
the  NEBRFD  specified.  At  that  small  angular  distance  from  the  specular  direction, 
however,  such  sensitivity  is  not  necessary.  It  is  necessary  at  the  other  extreme, 
at  the  maximum  baekscatter  angle.  Thus,  for  the  first  four  detectors,  we  have  used 
only  apertures,  but  for  the  next  six  we  have  used  lenses  to  collect  the  radiation 
to  obtain  the  sensitivity  specified. 

System  details 

A  schematic  diagram  of  the  system  is  shown  in  Fig.  4  A  photograph  is  shown  in 
Fig.  5  The  laser  shown  on  top  is  a  Hughes  Series  3000  lmW  helium-neon  laser  that 
has  its  output  at  0.6328um.  This  shines  on  the  folding  flat  mirror  to  the  right  of 
it  through  a  synchronous  chopper.  The  chopper  is  a  Bentham  218  Variable  Frequency 
Optical  Chopper.  The  beam  has  a  divergence  of  1.3mr  so  that  at  the  sample  shown 
below,  the  spot  has  spread  by  1.3mm,  and  is  therefore  14mm  in  diameter.  The  folding 
flat  is  of  our  own  manufacture.  The  specular  beam  enters  the  cylindrical  tube 
shown  where  there  are  two  crosshairs  that  constitute  an  aiming  device.  This  device 
allows  the  orientation  of  the  laser  and  the  detectors  to  the  local  surface  normal. 
It  establishes  the  angular  coordinate  system  to  the  extent  that  all  the  detectors 
are  now  in  the  plane  of  incidence  and  they  are  at  the  intended  angular  positions. 
The  sight  was  the  brainchild  of  A.  G.  DeBell  and  manufactured  by  us.  We  note  in 
passing  that  this  is  not  used  as  is  a  normal  rifle  sight.  Rather,  the  user  looks 
into  the  sight  first  from  the  front  and  then  from  the  rear,  each  time  at  an  angle 
to  the  optical  axis.  By  trial  and  error  he  insures  that  the  laser  beam  is  centered 
on  both  crosshairs.  This  way  he  accomplishes  the  alignment  without  losing  the  use 
of  his  sighting  eye. 

The  ten  silicon  detectors  are  all  model  UDT-455  Photops  detectors.  They  come 
complete  with  integral  preamplifiers.  They  were  mounted  in  structures  designed  and 
built  by  us  that  can  be  adjusted  to  view  the  spot  directly.  This  requires 
primarily  one  degree  of  freedom.  There  is  some  alignment  allowed  in  the 
out-of-plane  direction  so  that  we  can  be  sure  they  are  in  the  plane.  Figure  6  shows 
this  arrangement. 

Figure  7  is  a  schematic  of  the  preamplifier  that  is  integral  with  the  Photops 
detector.  These  detector  outputs  are  sampled  with  a  combination  of  hardware  and 
software.  The  hardware  consists  of  several  CMOS  chips:  2-IHSlOS's.  2-74HC244's. 
1-74HC373  ,  1-74HC74  ,  1-74HC138.  1-PGA200.  1-ICL7109  and  miscellaneous  components. 

The  master  controller  is  a  Radio  Shack  TRS-80  Model  100  used  for  its  properties  as 
a  terminal.  The  computer  serves  as  the  lockin  amplifier,  using  the  reference 
signal  from  the  chopper. 

Initial  checkout 


We  started  by  simply  checking  the  outputs  ag 
detectors.  We  used  Macor  as  a  represents 
standard,  and  we  used  a  mirror  available  in  the 
Early  tests  showed  that  detector  #1  which 
degrees  from  specular  gave  about  IV  signal  from 
about  2mV  signal:  the  noise  was  about  O.lmV  r 
as  being  white  noise.  A  detector  which  has  no 
specular  gave  respectively  20mV  and  5mV.  We 
data.  The  noise  level  is  comparable  to  a  BRDF 
lens  and  IE-2  without.  Some  tuning  should  lmpro 


ainst  two  samples  fo 
tive  of  a  highly  re 
lab  as  a  prototype 
has  a  lens  and  is 
the  Macor ,  while  t 
ms  and  showed  up  on 
lens  and  is  about  1 
can  draw  some  conclu 
of  about  IE-4  for  a 
ve  these  values  . 


r  two  different 
flective  diffuse 
specular  sample, 
located  about  80 
he  mirror  gave 
the  oscilloscope 
5  degrees  from 
sions  from  these 
detector  with  a 


Ca  l_i  b  rat  ion  tec  hn_i  qu  e 

The  system  needs  careful  handling  by  way  of  calibration.  Two  steps  in  general 
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will  be  taken.  The  first  is  to  use  a  diffuse  isotropic  sample  like  Macor.  It  will 
be  measured  by  our  laboratory  sea t t erome t er .  Then  the  channel  gains  will  be  set  in 
the  software  to  provide  the  same  BRDF  profile.  In  this  sense  then  the  instrument 
is  calibrated  with  Macor.  or  Halon  or  equivalent  as  a  reference.  Then  a  mirror 
sample  will  be  used  in  the  same  way,  first  on  the  laboratory  Instrument  and  then 
ith  the  portable  sea 1 1 eromet er  .  The  measurements  should  agree.  Future  use  will 
e  predicated  on  reca 1 1  bra t ion  against  the  isotropic,  diffuse  standard. 

Plans  for  the  future 


The  system  still  does  not  do  all  that  we  would  like  it  to  do. 
limited  in  its  capability  for  making  measurements  close  to  the  specul 
One  solution  to  this  problem  is  the  use  of  an  array  of  silicon  detect 
a  single  first  detector.  Of  course  this  requires  careful  cons 
calibration  of  the  beam  profile  in  the  area  of  the  specular  direct 
area  of  investigation  is  the  design  and  construction  of  a  similar  mod 
region  of  the  spectrum.  We  are  also  addressing  this  question. 
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APPENDIX  C 

Calibration  and  Operation  of  the  Portable  Scatterometer 

ASSEMBLY 


1.  Attach  the  crossrail  to  the  optics  rail  with  the  gusset  plate  and  install  the  two 
brackets  under  the  gusset  plate  on  each  side  of  the  optics  rail  (Fig.  1). 

2.  Attach  the  three  swivel  mounts  to  the  corners  of  the  optics  rail  and  crossrail  (Fig. 
2).  The  legs  of  the  tripods  should  be  retracted,  in  their  minimum  extension  position. 
Attach  the  tripods  to  the  swivel  mounts  by  inserting  the  threaded  stud  at  the  top  of 
the  tripod  into  a  matching  thread  on  the  bottom  of  the  swivel  mount  and  turning  the 
black  knurled  thumb  wheel  until  the  stud  is  firmly  secured  in  the  threaded  hole. 

3.  The  power  to  the  main  electronics  cabinet  should  be  off  during  assembly.  Plug 
the  laser  and  chopper  power  supplies  into  the  115  VAC  outlets  on  the  front  panel  on 
the  main  electronics  cabinet.  Connect  the  detector  power  cables  and  signal  cables  to 
the  front  panel  also  (Fig.  3);  each  cable  is  tagged  with  the  appropriate  detector 
number.  The  array  power  and  signal  cables  should  be  plugged  into  the  tenth  set  of 
sockets. 

4.  Connect  the  interface  cable  coming  from  the  bottom  of  the  TRS-80  model  100 
computer  to  the  interface  connector  on  the  rear  panel  of  the  main  electronics  cabi¬ 
net. 

CAUTION:  When  the  interface  cable  is  plugged  into  the  main  electronics  cabinet, 
care  must  be  taken  to  ensure  that  computer  is  never  powered  while  the  main  elec¬ 
tronics  cabinet  is  not.  FAILURE  TO  DO  SO  MAY  RESULT  IN  LOSS  OF  ALL 
PROGRAMS  AND  DATA  STORED  IN  RAM.  Therefore,  when  the  interface  cable 
is  plugged  in,  the  6  VDC  power  supply  for  the  computer  must  be  plugged  into  the 
115  VAC  time-delayed  outlet  indicated  in  Fig.  3.  Note  that  whenever  the  main 
electronics  cabinet  power  is  turned  off,  though,  the  computer  power  must  be  reset. 
To  reset  it,  turn  the  computer  power  switch  located  on  the  right  side  of  the  com¬ 
puter  off  and  then  on  again.  Once  reset  the  computer  will  be  powered  whenever  the 
main  electronics  cabinet  is  turned  on  again. 

5.  Connect  one  end  of  the  chopper  reference  cable  to  the  reference  output  connector 
on  the  chopper  power  supply  and  the  other  to  the  chopper  reference  connector  on 
the  front  panel  of  the  main  electronics  cabinet.  Connect  the  chopping  head  cable  to 
the  chopper  power  supply  also.  The  frequency  dial  of  the  chopper  power  supply 
should  be  set  at  about  28  for  a  5  Hz  chopper  frequency. 

6.  If  necessary,  loosen  the  laser  mount  screws  (Fig.  4)  and  remove  any  packing  mat¬ 
erial.  Position  the  laser  so  that  it  is  approximately  centered  over  the  two  mounts. 
The  laser  rotational  position  should  be  set  so  that  the  power  cord  is  emerging  from 
the  laser  body  at  the  bottom.  Re-tighten  the  laser  mount  screws. 
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OPERATION 


Alignment 

The  instrument  has  been  designed  to  measure  the  Bidirectional  Reflectance  Distri¬ 
bution  Function  (BRDF)  for  an  incidence  angle  of  45  degrees  (Fig.  4).  The  position 
of  the  laser,  turning  mirror,  and  linear  array  have  been  adjusted  so  that  when  the 
specular  beam  hits  the  center  of  a  removable  cross  hair  placed  in  the  slit  at  the  front 
of  the  array  sighting  tube  and  the  center  of  the  end  array  element  at  the  same  time, 
the  angle  of  incidence  is  45  degrees  and  all  detectors  are  in  the  plane  of  incidence. 
During  shipping  and  handling  of  the  instrument  the  position  of  these  items  may  have 
shifted  slightly.  The  alignment  can  be  checked  as  described  below  and  if  found  to 
be  out  of  alignment  perform  the  calibration  of  the  optical  system  described  in  the 
last  section  before  making  BRDF  measurements. 

1.  Loosen  the  tripod  post  lock  screws  at  the  top  of  each  tripod  (Fig.  2)  and  crank  up 
each  tripod  post  until  the  top  of  the  optics  rail  and  crossrail  is  about  25  inches  above 
the  floor  at  each  tripod.  Temporarily  tighten  all  tripod  post  lock  screws. 

2.  The  calibration  of  the  optical  system  is  based  on  the  center  being  defined  at  35.95 
inches  from  the  end  of  the  optics  rail  (Fig.  5).  Readjust  the  tripods  until  the  top  of 
the  optics  rail  at  the  center  is  25  inches  above  the  floor  on  both  sides  (points  A  and 
B  in  Fig.  5)  and  the  optics  rail  is  approximately  level.  Due  to  a  slight  twist  in  the 
rail  the  ends  may  be  slightly  higher  than  the  center.  A  circular  level  can  be  used  to 
get  the  rail  approximately  level  across  the  long  dimension  and  level  across  the  short 
dimension  everywhere  but  near  the  crossrail  (Fig.  5);  tripod  2  may  be  need  to  be 
slightly  higher  than  the  other  two.  After  adjustment  all  tripod  post  lock  screws 
should  be  tight. 

3.  Turn  on  the  laser  power  by  flipping  the  power  switch  on  the  rear  panel  of  the 
main  electronics  cabinet  (up  position).  The  chopper  power  supply  can  be  turned  off 
at  this  time.  The  axis  of  the  laser  beam  should  be  parallel  to  the  optics  rail  (Fig.  5). 
This  alignment  can  be  checked  by  setting  a  12-inch  machinist  square  upright  on  the 
optics  rail  near  the  laser  output  as  shown  in  Fig.6;  the  edge  of  the  scale  should  be 
two  inches  from  the  side  of  the  optics  rail.  A  piece  of  white  paper  placed  in  front 
of  the  scale  edge  can  be  used  to  determine  whether  or  not  the  scale  edge  bisects  the 
laser  spot.  Using  the  above  procedure  check  the  alignment  of  the  laser  beam  near 
the  turning  mirror  also.  The  axis  of  the  laser  beam  should  be  two  inches  from  the 
side  of  the  optics  rail  near  both  the  laser  output  and  the  turning  mirror. 

4.  Place  the  pentaprism  and  prism  housing  under  the  optics  rail  on  the  floor  so  that 
the  incident  beam  is  approximately  centered  on  the  1/2  inch  diameter  hole  in  the 
metal  cover  plate  of  the  housing  and  the  exit  beam  is  reflected  upward  toward  the 
array.  Remove  the  beveled  nylon  cross  hair  at  the  front  of  the  array  sighting  tube  if 
necessary.  Center  the  beam  spot  on  the  end  array  element  by  adjusting  the  penta¬ 
prism  position  and  tilt  as  described  below  until  the  maximum  signal  amplitude  is 
observed. 

The  signal  can  be  measured  by  connecting  the  array  signal  cable  to  a  voltmeter 
or  an  oscilloscope,  and  setting  the  switch  on  the  array  housing  to  "S"  (select  "S"  for 
specular  samples  and  "D"  for  diffuse  samples).  The  chopper  and  computer  power 
should  be  on  at  this  time.  When  the  computer  is  first  turned  on  a  directory  of  files 
is  displayed  with  the  cursor  over  BASIC.  Hit  ENTER  and  type  the  following 
sequence: 


TRITOO 

-2 


Figure  5. 


Figure  6. 


OUT  119,0  <ENTER> 
OUT  118,75  <ENTER> 


The  first  line  sets  the  gain  of  the  programmable-gain-amplifier  to  one.  The  second 
line  selects  which  array  element  preamp  is  switched  to  the  signal  cable;  in  this  case 
the  end  array  element  is  selected. 

Misalignment  of  the  beam  spot  on  the  vertical  axis  of  the  end  array  element  is 
caused  by  an  incorrect  position  of  the  prism  along  an  axis  common  to  the  plane  of 
incidence  and  the  base  of  the  prism  (axis  1  of  Fig.  7).  Translate  the  prism  along  this 
axis  until  the  maximum  signal  amplitude  is  observed.  Note  that  the  pentaprism  pro¬ 
duces  an  exit  beam  that  is  always  90  degrees  from  the  incident  beam  regardless  of 
the  tilt  of  the  prism  about  an  axis  perpendicular  to  the  plane  of  incidence  (axis  2  of 
Fig.  7),  so  no  adjustment  of  the  prism  about  this  axis  is  necessary. 

The  plane  of  incidence,  however,  is  dependent  upon  the  tilt  of  the  prism  about 
axis  1  of  Fig.  7.  Any  tilt  of  the  prism  about  this  axis  will  cause  horizontal  misalign¬ 
ment  of  the  beam  on  the  end  array  element.  Adjust  the  fine  pitch  screw  on  the 
prism  housing  until  the  maximum  signal  amplitude  is  obtained. 

The  pentaprism  may  also  be  rotated  about  an  axis  perpendicular  to  the  base  (axis 
3  of  Fig  7).  If  the  prism  is  rotated  slightly  about  this  axis,  the  exit  beam  will  be 
parallel  to  the  incident  beam  but  translated  sideways  (Fig.  8).  To  avoid  this  the  pen¬ 
taprism  should  be  rotated  about  axis  3  so  that  the  beam  incident  on  the  prism  and 
the  beam  reflected  back  toward  the  turning  mirror  by  the  prism  are  aligned  in  the 
same  vertical  plane  (Fig.  9). 

When  the  pentaprism  is  adjusted  so  that  the  maximum  signal  amplitude  from  the 
end  array  element  is  observed  and  the  incident,  exit,  and  reflected  beams  are  in  the 
same  vertical  plane,  check  the  signal  from  the  next  array  element  from  specular  to 
be  sure  it  is  smaller  then  the  signal  at  the  end  element.  To  do  this  type: 

OUT  118,74  <ENTER> 

which  switches  the  next  array  element  to  the  array  signal  cable.  When  the  switch  on 
the  array  housing  is  set  on  "S"  the  preamp  gain  is  a  factor  of  ten  greater  for  this 
array  element  than  for  the  end  element,  so  take  this  into  account  when  comparing 
the  two  signals.  Reconnect  the  array  signal  cable  to  the  front  panel  of  the  main 
electronics  cabinet. 

5.  Install  the  beveled  nylon  cross  hair  on  the  pins  next  to  the  slit  at  the  front  of  the 
array  sighting  tube  so  that  the  beveled  surface  faces  outward  and  the  cross  hair 
covers  the  slit  near  the  bottom  end.  If  the  instrument  is  correctly  aligned  the  beam 
spot  will  be  centered  on  the  cross  hairs. 

6.  Detector  voltage  data  from  the  Macor  sample  collected  by  running  the  BRDF  pro¬ 
gram  stored  in  the  TRS-80  model  100  computer  can  be  used  to  determine  whether  or 
not  the  detectors  need  realignment.  Replace  the  pentaprism  with  the  rectangular 
Macor  sample  and  mask  assembly  so  that  the  beam  spot  is  approximately  centered  on 
the  Macor  and  run  the  BRDF  program  as  directed  in  the  next  section,  "Running  the 
BRDF  Program." 

7.  The  detector  voltages  read  by  the  computer  with  the  Macor  sample  should  be 
within  10%  of  those  listed  in  Table  1. 


20 


M\c. 

C, 

l 


Table  1. 


Detector 

1 

2 

3 

4 

5 

6 

7 

8 
9 


Voltage  (volts) 
0.850 
1.35 
1.52 
1.56 
1.25 
0.985 
0.032 
0.029 
0.027 


If  they  are,  the  alignment  check  is  completed  and  the  scatterometer  is  ready  to  set 
up  for  measuring  the  BRDF  of  a  particular  mirror. 


Running  the  BRDF  Program 

When  the  computer  is  first  turned  on  a  directory  of  BASIC  programs,  text  files, 
and  data  files  will  be  displayed.  A  display  angle  control  located  on  the  right  side  of 
the  computer  can  be  adjusted  for  optimum  viewing  of  the  computer  display.  To 
return  to  the  directory  from  a  BASIC  program  or  text  file  hit  the  "F8”  key.  To  stop 
execution  of  a  program  at  any  time  press  the  shift  and  break  keys  simultaneously. 

Make  sure  a  text  file  called  KFACTR.DO  is  listed  in  the  directory  before  running 
the  BRDF  program.  To  run  the  program,  place  the  black  cursor  over  the  directory 
entry  PBRDF.BA  and  press  ENTER.  The  cursor  is  moved  by  pressing  one  of  the 
arrow  keys  in  the  upper  right  hand  corner  of  the  keyboard.  Type  in  the  sample  ID, 
date,  and  number  of  chopper  cycles  to  sample  for  each  detector  as  directed  by  the 
program.  Set  the  switch  on  the  array  housing  to  "S"  for  specular  samples  and  "D"  for 
diffuse  samples.  When  ENTER  is  pressed  after  entering  in  the  switch  setting,  the 
computer  will  read  the  detector  voltages  one  at  a  time  consecutively  starting  with  the 
one  farthest  from  specular.  It  is  recommended  that  all  room  lights  be  turned  off 
while  detector  signals  are  being  read.  The  computer  will  beep  once  after  each  det¬ 
ector  voltage  has  been  stored,  and  will  beep  twice  when  all  data  has  been  collected 
and  the  room  lights  can  be  turned  on  again.  If  desired,  save  the  data  in  a  data  file 
by  typing  in  a  file  name  and  whether  or  not  it  is  to  be  appended  to  an  existing  file 
as  directed  by  the  program  (appending  files  saves  directory  space).  The  data  can 
also  be  printed  out  on  a  line  printer. 


During  a  data  run  the  BRDF  program  may  display  one  of  two  error  messages. 
An  error  message  is  given  if  the  laser  output  fluctuates  more  than  5%  during  a  data 
run  or  when  the  program  autoranging  routine  has  set  the  amplifier  gain  to  one  (the 
lowest  gain  available)  and  the  signal  is  still  too  high.  Take  appropriate  steps  to  cor¬ 
rect  either  problem  before  running  the  program  again. 


BRDF  Measurements 


To  measure  the  BRDF  of  a  particular  mirror  it  is  necessary  to  align  the  instru¬ 
ment  in  reference  to  the  surface  being  measured.  The  turning  mirror,  laser,  array, 
and  detectors  should  not  be  readjusted  at  this  point;  adjust  only  the  tripods  as  out¬ 
lined  below. 

1.  Place  the  instrument  over  the  mirror  to  be  measured  as  shown  in  Fig.  4.  Loosen 
the  tripod  post  lock  screws  and  adjust  the  tripods  so  that  the  top  of  the  optics  rail  at 
the  center  is  about  24  inches  above  the  mirror  surface  and  approximately  parallel  to 
it,  and  the  specular  beam  is  incident  on  the  slit  at  the  front  of  the  array  sighting 
tube.  Temporarily  tighten  all  tripod  post  lock  screws. 

2.  Remove  the  beveled  nylon  cross  hair  at  the  front  of  the  array  sighting  tube  if 
necessary.  Center  the  beam  on  the  end  array  element  by  adjusting  the  tripods  until 
the  maximum  signal  amplitude  is  obtained.  The  procedure  for  measuring  the  signal 
voltage  is  described  in  step  four  of  the  "Alignment"  section  above. 

Alignment  along  the  vertical  axis  of  the  erd  array  element  is  done  by  adjusting 
tripod  1  (see  Fig.  5),  and  along  the  horizontal  axis  by  adjusting  tripods  2  and  3. 
Note  if  the  top  of  the  optics  rail  is  not  the  correct  height  above  the  sample  surface, 
it  may  be  necessary  to  adjust  all  three  tripods  by  the  same  amount  to  get  the  right 
height  and  still  keep  the  optics  rail  parallel  to  thf  sample  surface. 

3.  When  the  beam  spot  is  centered  on  the  end  array  element,  install  the  beveled 
nylon  cross  hair  so  that  the  beveled  surface  faces  outward  and  the  cross  hair  covers 
the  slit  near  the  bottom  end.  If  the  specular  beam  is  not  centered  on  the  cross  hair, 
readjust  the  tripods.  Remove  the  cross  hair  and  recheck  the  alignment  on  the  end 
array  element.  Continue  in  this  manner  until  the  specular  beam  is  centered  on  both 
the  cross  hair  and  the  end  array  element.  After  adjustment  all  tripod  post  lock 
screws  should  be  tight. 

4.  Run  the  BRDF  program  as  described  under  "Running  the  BRDF  Program". 
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CALIBRATION  OF  THE  OPTICAL  SYSTEM 


Calibration  of  the  optical  system  includes  setting  up  the  scatterometer  for  meas¬ 
urements  at  a  45-degree  incidence  angle,  aligning  the  axis  of  the  linear  array  in  the 
plane  of  incidence,  aligning  the  detectors  for  maximum  output,  and  determining 
BRDF  gain  factors.  The  gain  factors  for  each  detector  and  array  element  are  found 
by  dividing  BRDF  values  obtained  from  Macor  on  a  laboratory  scatterometer  of 
known  performance  by  the  voltages  measured  from  the  Macor  on  the  portable  scat¬ 
terometer.  The  portable  scatterometer  BRDF  program  calculates  the  BRDF  of  a 
given  surface  by  multiplying  these  gain  factors  by  the  measured  signals  from  that 
surface.  The  detector  angles  were  calculated  by  assuming  the  top  of  the  optics  rail 
was  24  inches  above  the  sample  surface  and  defining  the  center  of  the  optics  rail  at 
35.95  inches  from  the  end  of  the  rail  (Fig.  10). 

1.  All  of  the  surfaces  used  in  the  calibration  procedure  are  mounted  so  that  the  sur¬ 
face  is  one  inch  above  the  floor,  i.e.,  the  one  inch  target  block,  the  pentaprism,  and 
the  Macor.  Therefore,  it  is  necessary  to  set  the  top  of  the  optics  rail  25  inches 
above  the  floor  so  that  it  is  24  inches  above  the  sample  surface.  Adjust  the  tripods  as 
described  in  steps  one  and  two  in  the  "Alignment"  section  above. 

2.  Check  the  alignment  of  the  laser  as  directed  in  step  three  in  the  "Alignment"  sec¬ 
tion  above.  If  the  laser  is  out  of  alignment  loosen  the  laser  support  block  screws 
located  on  the  bottom  of  the  optics  rail.  Re- position  the  support  blocks  until  the 
beam  is  parallel  to  the  optics  rail,  and  re-tighten  the  support  block  screws. 

3.  For  a  45-degree  angle  of  incidence  at  the  sample  surface  the  turning  mirror  must 

be  positioned  so  that  the  distance  from  the  center  to  the  point  where  the  laser  hits 
the  turning  mirror  is  29.25  inches  as  shown  in  Fig.  10.  Remove  the  turning  mirror 
housing  from  the  optics  rail  for  easier  access  to  the  turning  mirror.  The  turning 
mirror  base  is  adjustable  along  the  axis  of  the  optics  rail  and  is  held  in  place  by  two 

screws  in  the  slotted  base  plate.  Position  the  base  so  that  the  outer  edge  is  5.15 

inches  from  the  end  (Fig.  10)  and  the  base  is  approximately  parallel  to  the  rail,  and 

tighten  the  base  plate  screws. 

4.  The  laser  and  chopper  power  can  be  turned  off  at  this  time.  Position  the  plumb 
bob  bracket  so  that  its  center  is  at  the  center  of  the  optics  rail  on  both  sides  (points 
A  and  B  in  Fig.  5).  Place  the  one  inch  thick  target  block  on  the  floor  so  that  the 
point  of  the  plumb  bob  is  centered  on  the  cross  hair.  Remove  the  plumb  bob  once 
the  block  is  positioned. 

5.  Turn  the  laser  power  on  and  adjust  the  turning  mirror  until  the  beam  is  centered 
on  the  cross  hair  on  the  one  inch  target  block.  The  mirror  is  adjustable  by  means  of 
fine  pitch  adjustment  screws  on  the  rear  of  the  mirror  mount  (Fig.  11).  The  beam 
should  pass  through  the  top  of  the  optics  rail  at  24  inches  from  the  center  (Fig.  10) 
and  2  inches  from  the  side;  if  not  some  slight  adjustments  of  the  turning  mirror  may 
be  necessary.  Install  the  turning  mirror  housing. 

6.  Replace  the  one  inch  target  block  with  the  pentaprism  and  prism  housing  so  that 
the  incident  beam  is  approximately  centered  on  the  1/2- inch  diameter  hole  in  the 
metal  cover  plate  of  the  housing,  and  the  exit  beam  is  reflected  upward  toward  the 
array.  Adjust  the  position  and  tilt  of  the  prism,  as  discussed  in  step  four  of  the 
"Alignment"  section  above,  until  the  exit  beam  passes  through  the  top  of  the  optics 
rail  24  inches  from  the  center  at  the  array  end  of  the  rail  and  2  inches  from  the 
side. 


Figure  10. 


7.  Install  the  beveled  nylon  cross  hair  on  the  pins  next  to  the  slit  at  the  front  of  the 
array  sighting  tube  so  that  the  beveled  surface  faces  outward  and  the  cross  hair 
covers  the  slit  near  the  bottom  end.  The  rotational  position  of  the  array  assembly 
within  the  plane  of  incidence  is  held  in  place  by  two  nuts  at  the  pivot.  Adjust  the 
array  assembly  so  that  the  sighting  tube  is  about  45  degrees  from  the  horizontal  and 
temporarily  tighten  the  pivot  nuts. 

8.  The  array  assembly  is  held  in  place  on  the  optics  rail  by  two  screws  in  the  slotted 
base  plate.  Due  to  the  twist  in  the  optics  rail  it  may  be  necessary  to  put  shims  under 
the  base  plate  so  that  the  array  assembly  will  be  level.  If  it  is  not  level,  the  axis  of 
the  array  will  not  be  parallel  to  the  plane  of  incidence  which  is  perpendicular  to 
level  surfaces.  A  circular  level  placed  on  the  optics  rail  may  be  used  to  determine 
how  much  compensation  is  needed  to  make  the  array  assembly  level.  Position  the 
base  plate  with  the  shims  under  it  so  that  the  beam  is  approximately  centered  on  the 
cross  hair  and  temporarily  tighten  the  base  plate  screws. 

9.  Remove  the  beveled  nylon  cross  hair  and  measure  the  signal  amplitude  of  the  end 
array  element  using  the  procedure  described  in  step  four  of  the  "Alignment"  section 
above  (the  chopper  should  be  turned  on  at  this  time).  A  lock-in  amplifier  with  the 
reference  signal  obtained  from  the  chopper  can  be  used  instead  of  the  oscilloscope  or 
voltmeter.  Slightly  loosen  the  pivot  nuts  of  the  array  assembly,  adjust  the  rotational 
position  until  the  maximum  signal  from  the  end  array  element  is  observed,  and  re¬ 
tighten  the  pivot  nuts.  Slightly  loosen  the  base  plate  screws,  rotate  the  array  assem¬ 
bly  side  to  side  for  maximum  signal,  and  re-tighten  the  base  plate  screws. 

10.  Replace  the  beveled  nylon  cross  hair.  If  the  beam  is  not  centered  on  the  cross 
hair  slightly  loosen  the  base  plate  screws  and  re-position  the  array  assembly. 
Remove  the  cross  hair  and  repeat  step  nine.  Continue  in  this  manner  until  the  beam 
is  centered  on  both  the  cross  hair  and  the  end  array  element.  Then  check  the  signal 
from  the  next  array  element  to  be  sure  it  is  smaller  than  the  signal  from  the  end  ele¬ 
ment  as  discussed  in  step  four  of  the  "Alignment"  section  above. 

11.  The  chopper  frequency  should  be  5  Hz.  If  not,  adjust  the  frequency  dial  on  the 
chopper  power  supply;  a  setting  of  approximately  28  produces  a  5  Hz  chopper  fre¬ 
quency.  Reconnect  the  array  signal  cable  to  the  front  panel  of  the  main  electronics 
cabinet. 

12.  Adjust  the  nine  detector  assemblies  for  maximum  signal  amplitude  at  each  output 
as  follows.  For  each  detector,  disconnect  the  signal  cable  from  the  electronics  cabi¬ 
net  and  connect  it  to  an  oscilloscope  or  lock-in  amplifier  (the  chopper  should  be  on 
at  this  time).  Loosen  the  four  screws  holding  the  detector  in  place  on  the  optics  rail, 
adjust  the  position  of  the  detector  assembly  until  the  maximum  signal  is  observed, 
and  re-tighten  the  screws.  Loosen  the  pivot  nuts  and  screws  slightly  (Fig.  12)  so  that 
the  detector  housing  is  tumable  but  not  loose.  Slowly  rotate  the  detector  housing  to 
obtain  the  maximum  signal,  and  re-tighten  the  pivot  nuts  and  screws.  Reconnect  the 
detector  signal  cable  to  the  front  panel  of  the  main  electronic  cabinet. 

13.  When  the  signal  of  all  nine  detectors  has  been  maximized,  run  the  BRDF  pro¬ 
gram  as  described  above  with  the  Macor  sample  and  save  the  data.  Return  to  the 
directory  and  place  the  cursor  over  KGENER.BA  and  hit  ENTER.  This  program 
calculates  the  gain  factors  for  each  detector  and  array  element  and  stores  them  in  the 
KFACTR.DO  file  used  by  the  portable  scatterometer  BRDF  program.  It  requires  the 
voltages  measured  by  each  detector  and  array  element  from  the  Macor  sample  and 
these  may  be  enter  from  the  newly  created  data  file  or  manually  from  the  keyboard. 
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Figure  12. 
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14.  Rerun  the  BRDF  program  with  the  Macor  sample  with  the  new  gain  factors. 
The  BRDF  should  be  about  I/*  at  all  detector  and  array  element  angles. 
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Tucson,  Arizona 

Abstract 


A  helium-neon  portable  scatterometer  that  vas  built  and  described 
previously  has  been  modified.  The  intent  vas  to  obtain  better 
data  at  angles  very  close  to  specular.  Ve  incorporated  a  tvelve- 
element,  linear  array  of  silicon  detectors  in  place  of  the  first 
detector  near  the  specular  reflection  direction,  and  removed  the 
collimator  to  reduce  the  beam  diameter.  Attention  vas  paid  to 
scattered  light  from  the  instrument.  A  nev  scheme  to  orient  the 
instrument  to  the  sample  surface  using  the  array  and  the  specular 
beam  is  described.  Design  details  of  the  detection  system  capa¬ 
ble  of  detecting  signals  ranging  over  five  orders  of  magnitude 
and  utilizing  lock-in  detection  software  are  given.  Results  shov¬ 
ing  that  the  modification  improved  the  scatterometer  performance 
so  that  a  Bidirectional  Reflectance  Distribution  Function  (BRDF) 
value  of  1  sr~  at  one  degree  from  specular  can  be  measured  with 
reasonable  accuracy  for  a  specular  sample  are  discussed. 


Introduction 


As  optical  systems  become  more  sophisticated  the  effects  of  light 
scatter  from  optical  components  on  system  performance  becomes 
Increasingly  Important,  especially  at  angles  very  close  to  specu¬ 
lar.  A  helium-neon  portable  scatterometer  capable  of  measuring 
the  Bidirectional  Reflectance  Distribution  Function  (BRDF),  that 
character lzes  scatter  from  component  surfaces,  vas  built  and 
described  in  an  earlier  paper  .  The  previous  configuration  is 
shown  in  Fig.  1  .  The  output  of  a  Hughes  Series  3000  helium-neon 
laser  (0.6328aub)  vas  collimated,  chopped  at  a  frequency  of  5  hz, 
and  reflected  to  the  sample  surface  by  a  folding  flat  shown  on 
the  right.  Instrument  orientation  to  the  sample  surface  was  done 
by  centering  the  specular  beam  in  the  sight  at  the  far  left.  Ten 
hybrid  silicon  photodetector /ampl if iers  each  mounted  at  the  end 
of  sighting  tubes  were  placed  at  various  angles  between  the 
source  and  specular  beam;  the  first  four  near  specular  were 
equipped  with  apertures,  the  other  six  vith  lenses.  A  Radio  Shack 
TRS-80  Model  100  Portable  Computer  was  used  to  gather  and  process 
the  signals  from  the  ten  detectors.  This  portable  instrumenta¬ 
tion  was  capable  of  measuring  scatter  from  diffuse  and  specular 
samples  up  to  two  meters  in  diameter  at  0.6328pm  for  a  43  degree 
incident  angle  over  a  fairly  large  range  of  scatter  angles.  How¬ 
ever,  measurements  were  limited  to  a  minimum  angle  from  specular 
of  1.5  degrees  which  corresponds  to  a  fi-fi  of  0.019  where  i® 
the  surface  spatial  frequency  times  the  wavelength  over  2  . 

This  paper  discusses  the  addition  of  a  twelve-element,  linear 
array  of  silicon  detectors  placed  near  the  specular  beam  to 
obtain  better  BRDF  data  in  the  region  near  specular. 

Modification  Details 


Ve  modified  the  scatterometer  by  replacing  the  sight  and  the 
detector  closest  to  specular  with  the  array  assembly.  The  modi¬ 
fication  is  shown  in  Fig.  2.  The  collimator  at  the  output  of  the 
laser  was  also  removed  to  reduce  the  specular  beam  diameter  at 
the  array  from  14mm  to  4.5mm.  Design  considerations  that  arise 
as  a  result  of  the  modification  are  the  field  of  view  of  each 
array  element,  the  scattered  light  from  the  instrument,  instru¬ 
ment  alignment  to  the  sample  surface  now  that  the  sight  has  been 
removed,  and  the  need  for  a  fairly  large  detector  dynamic  range. 

The  field  of  view  of  the  array  elements  is  defined  by  a  27mm  X 
3 . 6mm  slit  placed  on  the  end  of  the  array  sighting  tube.  This 
provides  a  rectangular  field  of  view  at  the  sample  with  a  minimum 
dimension  of  10.4mm  in  the  direction  normal  to  the  plane  of  inci- 
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FIGURE  1  Previous  scatterometer  configuration 


dence  for  each  array  element.  This  field  of  view  is  fairly  large 
but  the  elgnal-to-nolee  ratios  near  specular  should  generally  be 
high,  and  attention  was  paid  to  scattered  light  from  the  instru¬ 
ment.  To  minimize  this  scatter  we  replaced  the  folding  flat  with 
a  smoother  mirror  and  covered  all  reflective  parts  of  the  instru¬ 
ment. 

Since  the  detector  array  replaced  the  sight  UBed  for  instrument 
orientation  to  the  sample  surface,  we  incorporated  the  array  in 
the  new  alignment  scheme.  Previously,  the  specular  beam  was  cen¬ 
tered  on  crosshairs  at  both  ends  of  a  eight  that  was  set  at  45 
degrees  in  the  instrument  angular  coordinate  system.  This  alignr 
ment  insured  that  the  angle  of  incidence  on  the  sample  was  45 
degrees  and  the  detectors  were  in  the  plane  of  incidence  at  spe¬ 
cified  angles.  The  new  arrangement  is  based  on  the  same  prin¬ 
ciple.  The  array  assembly  is  set  at  45  degrees  and  the  beam  is 
centered  on  a  removable  crosshair  placed  in  the  slit  at  the 
front.  At  the  other  end  of  the  sighting  tube  the  beam  is  cen¬ 
tered  on  the  end  array  element. 

Detection  System 

The  detector  arra;  '-..rdware  includes  a  preamplifier  for  each  of 
the  twelve  arr^y  elements,  an  Analog  Devices  7506  1 6-channel 
analog  multiplex'..  ,  a  Burr  Brown  programmable-gain-amplif  ier ,  and 
an  Intersil  "109  integrating  A/D  converter.  A  block  diagram  is 
shown  in  Pig.  3.  Early  tests  showed  signals  gathered  from  the 
array  with  a  highly  diffuse  sample  and  with  a  low  scatter  mirror 
ranged  over  five  orders  of  magnitude.  This  range  was  too  large 
to  maintain  both  linearity  over  the  entire  range  and  accuracy 
from  the  A/D  converter  at  low  signal  levels.  Ve  designed  the 
preamplifiers  to  maintain  linear  operation  over  the  desired  range 
and  used  the  programmable-gain-amplifier  to  provide  more  gain  at 
low  levels.  In  addition,  a  switch  is  provided  to  manually 
increase  the  gain  on  the  pr eampl if iers  for  the  first  three  array 
elements  off  specular  when  measuring  diffuse  samples  and  decrease 
the  gain  for  specular  samples.  Control  of  the  A/D  converter,  the 
multiplexer,  and  the  programmable-gain-amplifier  is  done  via 
ports  available  in  the  TRS-80  Model  100  computer. 

Ve  used  a  software  approach  to  lock-in  detection  to  minimize  cost 
and  increase  instrument  portability.  Each  detector  signal  is 
sampled  over  several  chopper  cycles.  The  computer  monitors  the 
chopper  and  repetitively  adds  signals  sampled  while  the  chopper 
is  open  and  subtracts  while  the  chopper  is  closed.  The  software 
ensures  that  the  number  of  additions  equals  the  number  of  sub¬ 
tractions  and  thus  the  noise  averages  to  zero  while  the  signal 
just  keeps  adding.  The  signal  total  is  then  divided  by  the  num- 
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FIGURE  3  Block  diagram  of  the  detection  system 


ber  of  cycles  sampled.  The  lock-in  portion  of  the  system  soft¬ 
ware  is  written  in  3085  assembly  language  in  order  to  convert  and 
add  or  subtract  signals  within  a  half  cycle  of  the  chopper 
frequency. 


Results 


Demonstration  of  performance  at  near-specular  angles  and  rela¬ 
tively  low  BRDF  values  is  difficult.  Ve  approached  this  in  sev¬ 
eral  steps.  The  first  was  to  measure  Macor ,  essentially  a  per¬ 
fect  isotropic  sample.  The  result  as  shovn  in  Fig.  4  is  1/  "all 

across  the  board."  The  data  points  between  P-P  *  0.001  and  P-P 

0  0 

*  0.015  were  obtained  from  the  array  while  the  other  nine  points 
were  from  the  remaining  detectors.  The  next  step  was  to  measure 
one  of  our  best  mirrors,  a  mirror  made  some  time  ago  by  Jaeger, 
and  compare  it  with  the  instrument  profile.  The  profile  was 
measured  by  setting  the  instrument  in  a  transmission  arrangement 
with  no  sample  as  shown  in  Fig.  5.,  and  is  represented  by  the 
crosses  in  Fig.  6.  There  is  a  strong  indication  that  the  profile 
drops  rapidly  near  P-P  *  0.01.  The  points  further  out  are  the 
noise  equivalent  BRDF  (&EBRDF);  the  last  six  points  were  obtained 
from  detectors  equipped  with  lenses  producing  a  larger  collecting 
area  and  thus  a  lower  NEBRDF.  Under  the  conditions  of  measure¬ 
ment  (45  degrees  incident,  0.6328^m)  the  mirror  still  looks  good 
as  shown  by  the  circles  in  Fig.  6.  Notice  that  from  P-P  *  0.045 
(3-5  degrees)  to  larger  angles  the  BRDF'  is  a  straight  line  with  a 
slope  of  about  -2.  This  agrees  with  our  expectations.  Starting 
at  P-P  =  0.015  (1.2  degrees)  the  curve  is  higher  and  no  longer 
fits  tfie  6ame  slope,  bv^t  closely  matches  the  instrument  profile. 
It  ha6  a  value  of  1  sr  at  one  degree.  Next,  a  second  inferior 
mirror  was  measured  and  compared  with  the  instrument  profile.  As 
shown  in  Fig.  7  this  curve  is  higher  with  a  slope  of  about  -2 
that  rolls  over  at  0.012  and  crosses  the  instrument  profile  at 
around  0.004  (0.3  degrees),  although  there  may  be  some  contribu¬ 
tion  starting  at  0.01  (0.8  degrees). 

Conclusion 


The  new  configuration  reduced  the  1 instrument  background  profile 
to  the  point  where  a  BRDF  of  Isr  can  be  measured  with  reason¬ 
able  accuracy  at  one  degree  from  specular. 
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FIGURE  5  Transmission  arrangement  for  measuring  the  instrument  profile 


FIGURE  7  Comparison  of  a  second  inferior  mirror  BRDF  with  the  instrument  profile 
(O-mirror,  X-profile) 
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A  10.6  fim  Portable  Scatterometer  for  Optical  Shop  Use 
Steven  J.  Wein  and  William  L.  Wolfe 
INTRODUCTION 

We  have  bread-boarded  and  tested  a  design  for  a  10.6-/im  portable 
scatterometer  which  can  test  mirrors  up  to  1.7  meters  in  diameter  with  the  emphasis 
on  obtaining  near-specular  measurements.  We  have  obtained  and  verified  a  design 
that  can  measure  scattering  from  2°  off  specular  to  80°  off  specular.  Near  specular 
the  NEBRDF  (noise  equivalent  bidirectional  reflectance  distribution  function)  is 
10~4  sr-1  and  far  off  specular  the  NEBRDF  is  approximately  10-6  sr-1. 

The  system  consists  of  a  7.5-W  C02  laser  whose  beam  is  expanded  and 
focused  by  a  zinc  selenide  beam  expander.  A  He-Ne  laser  beam  is  coalligned  with 
the  infrared  beam  for  aid  in  alignment  and  safety  purposes.  We  utilize  pyroelectric 
radiometers  to  detect  the  scattered  radiation. 

The  configuration  depicted  in  Figure  1  was  chosen  on  the  merits  of  portability 
and  compactness.  The  requirements  of  portability  place  special  constraints  on  the 
geometry  of  the  set-up  which  also  affect  the  near  specular  performance  (which  will 
be  subsequently  discussed).  In  addition,  this  instrument  is  intended  to  be  able  to  test 
many  different  mirrors  with  a  wide  range  of  radii  of  curvature.  Therefore,  in 
general,  we  cannot  rely  upon  the  power  of  the  test  surface  to  aid  in  the  reduction  of 
the  focal  spot  diameter.  We  are  able  to  make  BRDF  measurements  to  2°  from 
specular  when  testing  a  flat  mirror.  When  testing  a  mirror  with  positive  power,  we 
can  make  measurements  inside  2°  depending  on  the  power  of  test  surface.  On  the 
other  hand,  we  cannot  quite  obtain  2°  from  specular  when  testing  a  surface  with 
negative  power. 


OPTICAL  PATH  LENGTH 

The  above  design  results  in  two  conflicting  requirements  on  the  length  of 
section  B.  The  first  requirement  is  to  minimize  the  focal  spot  diameter.  In  order  to 
do  this,  the  system  f /#  must  be  reduced;  hence,  section  B  must  be  reduced  in  length. 
The  second  requirement  is  to  minimize  the  system  scattering.  To  do  this  the 
illuminated  portion  of  the  turn  mirror  must  be  outside  the  field  of  view  (FOV)  for 
all  the  radiometers.  Especially  for  the  near  specular  radiometer  this  requires  an 
increase  in  the  length  of  section  B.  On  paper,  there  is  a  region  that  satisfies  both  of 
these  conditions.  For  our  geometry,  both  conditions  are  theoretically  satisfied  for  1° 
off  specular  when  section  B  of  Figure  1  is  1.94  meters  or  slightly  greater. 
(Unfortunately,  this  is  rather  large  to  satisfy  the  portability  requirement.)  Both 
conditions  are  satisfied  for  1.8°  off  specular  when  section  B  is  approximately  60  cm. 
However,  this  assumes  that  the  baffle  of  the  radiometer  stops  all  off-axis  radiation 
and  that  the  specular  beam  falls  off  as  a  Gaussian  beam  over  orders  of  magnitude 
from  the  peak  irradiance.  Both  are  dubious  assumptions.  In  practice  we  have  found 
that  the  roll-off  of  the  beam  profile  is  sufficiently  low  and  the  baffle  of  the 
radiometer  is  sufficiently  "black"  at  2°  from  specular  to  obtain  BRDF  measurements 
of  10~4  sr-1. 

In  addition,  section  A  should  be  short  enough  to  yield  an  acceptable  signal-to- 
noise  ratio  (SNR)  for  a  BRDF  of  IO-4  sr-1.  However,  section  A  should  also  be  long 
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enough  that  the  separation  between  the  specular  beam  and  2°  is  large  enough  for  the 
specular  irradiance  to  be  below  the  scattered  irradiance. 

NEAR  SPECULAR  BAFFLE 

Off-axis  stray  radiation  inside  the  radiometer  is  perhaps  the  predominant 
problem.  Especially,  but  not  exclusively,  at  two  degrees  off  the  axis  of  the 
radiometer  is  a  strong  noise  contribution  due  to  the  system  scattering  and  the 
specular  beam. 

If  we  utilize  a  long  baffle  radiometer  we  can  forego  the  need  for  a  higher  D* 
detector.  To  obtain  a  sufficiently  small  FOV  for  the  radiometer  such  that  it  does  not 
see  the  "illuminated"  portion  of  the  turn  mirror  (which  we  are  defining  to  be  the 
1/e4  irradiance  points)  we  need  to  close  down  the  entrance  aperture  of  the  baffle. 
However,  for  short  baffles  the  aperture  becomes  sufficiently  small  such  that  the 
power  transfer  from  the  sample  to  the  detector  yields  a  signal  to  noise  less  than  one. 
The  solution  to  this  is  to  determine  the  appropriate  aperture  large  enough  to  allow  a 
sufficient  SNR  then  extend  the  baffle  length  to  reduce  the  FOV  to  be  sufficiently 
narrow.  A  3-ram  aperture  at  the  end  of  a  25.4-cm  baffle  meets  the  criteria 
(Figure  2).  We  placed  additional  fins  inside  the  baffle  as  indicated  above  such  that 
the  detector  cannot  directly  see  the  sides  of  the  baffle.  This  did  reduce  the  off-axis 
stray  radiation  rejection  but  not  sufficiently  to  make  measurements  at  one  degree  off 
specular. 
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Figure  2.  Baffle  schematic. 


BEAM  EXPANDER/FOCUSING  ELEMENT 

We  must  significantly  underfill  the  exit  pupil  of  the  system  to  maintain  the 
Gaussian  roll-off  of  the  focal  spot.  The  truncated  beam  will  very  nearly  resemble  a 
Gaussian  beam  until  the  point  of  truncation  where  we  obtain  ringing  in  the 
diffraction  pattern  that  falls  off  as  1/r3. 

The  exit  pupil  diameter  of  our  zinc  selenide  beam  expander  is  23  mm  and  the 
1/e2  diameter  of  the  Gaussian  beam  at  the  exit  pupil  is  7.5  mm.  Hence,  the  beam  is 
truncated  just  below  the  1/e18  level  and  theoretically  maintains  the  Gaussian  roll-off 
until  just  reaching  the  scattering  level.  This  was  verified  by  first  obtaining  a  signal 
voltage  due  to  a  gold  sandpaper  Lambertian  surface.  Then  an  rms  noise  voltage  was 
obtained  while  the  laser  was  operating  and  all  the  optical  and  electrical  noise 
mechanisms  where  active.  We  obtained  a  ratio  of  3000  between  the  two 
measurements  which  is  what  was  required  to  obtain  a  10-4  BRDF  measurement  at 
that  location. 

Spatial  filtering  of  the  beam  with  a  Keplerian  type  beam  expander  improves 
the  system  performance  such  that  we  are  able  to  make  measurements  to  1.60°  off 
specular.  However,  the  diameter  of  the  pinhole  is  very  important.  If  the  pinhole  is 
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too  small,  the  diffraction  from  the  pinhole  can  fill  the  exit  pupil.  If  the  exit  pupil  is 
fully  illuminated  with  a  plane  wave  for  example,  the  system  point  spread  function 
will  be  the  Airy  diffraction  pattern  which  does  not  roll  off  quickly  enough  to  make 
near-specular  measurements.  We  experimentally  determined  that  the  optimum  spatial 
filter  diameter  for  this  particular  application  is  approximately  3.5  times  the  diameter 
of  the  focal  spot  at  the  1/e2  points.  Spatial  filters  with  smaller  diameters  produced 

an  increase  in  specular  irradiance  in  the  region  from  one  to  two  degrees  from 

specular. 

DETECTION  SYSTEM 

This  system  utilizes  the  Molectron  PI -72  and  Pi -75  lithium  tantalate 

pyroelectric  detectors  with  built-in  JFET  preamplifiers.  These  are  low-noise 

pyroelectric  detectors  with  an  NEP  of  0.5  nW.  The  use  of  pyroelectric  detectors 
facilitates  portability  and  ease  of  operation  since  cooling  is  not  required.  The 
detectors  are  operated  in  a  voltage  gain  configuration  and  ac  coupled  into  a  low  noise 
Burr-Brown  OPA111  operational  amplifier  with  feedback  set  for  a  closed  loop  gain  of 
500.  This  yields  an  acceptable  signal  level  for  the  subsequent  A/D  conversion  and 
signal  processing. 


MECHANICAL  CHOPPING  SCHEME 

The  chopping  of  the  beam  by  the  chopper  blade  yields  a  temporal  convolution 
between  the  shape  of  the  slots  in  the  chopper  blade  and  the  beam  cross  section.  (See 
Figure  3.)  Thus,  we  can  modify  the  temporal  signal  as  seen  by  the  pyroelectric 
detectors  by  appropriately  choosing  the  width  and  shape  of  the  slots  in  the  chopper 


blade.  At  maximum  responsivity  our  pyroelectric  detectors  are  operating  in  the 
integrating  region  of  their  Bode  plot.  Utilizing  phase  sensitive  lock-in  detection  we 
would  like  as  much  as  possible  of  the  detector  output  to  be  in  the  fundamental 
harmonic.  Hence,  we  would  like  the  temporal  signal  to  be  .%s  close  to  a  sinusoid  as 
possible.  By  having  the  width  of  the  slots  (in  the  direction  normal  to  the 


measurement  plane)  approximately  equal  to  the  1/e2  diameter  of  the  Gaussian  beam, 
the  convolution  yields  a  temporal  signal  that  is  very  nearly  sinusoidal  when  viewed 
on  an  oscilloscope.  (The  other  dimension  of  the  slot  is  necessarily  large  to  avoid 
additional  unwanted  diffraction  in  the  measurement  plane.) 

All  results  described  above  were  obtained  using  the  Optical  Sciences  Center 
Infrared  Facility  Laboratory.  Thus,  we  have  been  able  to  show  with  laboratory 
components  an  instrument  configuration  that  has  an  instrumental  NEBRDF  of  10-4  sr"1 
at  2°  from  specular.  The  next  step  is  the  assembly  of  the  components  into  an 
operational  instrument. 
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Scattering  data  taken  at  0.6328.  3.39.  and  10.6  Mm  are  compared  with  theoretical  calculations  based  on  vector  and 
scalar  theories  of  scattering.  Methods  are  described  lor  obtaining  the  rms  surface  height  using  various  degrees  of 
dependence  on  the  models  of  scattering. 


INTRODUCTION 

Scattering  from  almost  smooth  surfaces  is  important  in  some 
optical  systems  and  in  certain  radar  problems.  These  sur¬ 
faces,  with  roughness  less  than  a  wavelength,  are  realized  by 
well-made  mirrors  in  optics  and  relatively  calm  ocean  surfaces 
at  longer  wavelengths.  Calculations  of  the  amount  and  dis¬ 
tribution  of  the  scattering  have  been  made  based  on  both  the 
scalar  and  the  vector  approaches  to  electromagnetic  boundary 
value  problems.  The  scalar  theory,  based  on  a  scalar  dif¬ 
fraction  integral,  has  been  developed  by  Davies.1  Beckmann,2 
and  Harvey.3  The  vector  theory,  based  largely  on  perturba¬ 
tion  or  variational  techniques  or  the  Stratton-Chu-Silver 
diffraction  integral,4-5  has  been  developed  by  many  authors.® 
The  perturbation  technique  was  originated  by  Rayleigh.7 
extended  by  Rice,  and  used  for  the  solution  of  both  optical  and 
radar  problems. 

Although  measurements  have  been  made  over  the  years  of 
scattering*  at  various  wavelengths9  and  for  surfaces  of  various 
roughnesses,  there  seems  to  have  been  no  systematic  com¬ 
parison  of  the  theory  and  the  data. 

We  have  made  measurements  of  the  angular  distribution 
of  scattering  for  three  different  angles  of  incidence,  five  dif¬ 
ferent  wavelengths,  and  several  different  polarization  ar¬ 
rangements.  all  for  «  single  sample  of  smooth  fused  silica  v,ith 
an  evaporated  aluminum  coating.  We  compared  these  data 
with  the  vector  and  the  scalar  models,  and  we  drew  some 
tentative  conclusions  about  the  agreement  of  theory  with 
measurements  based  on  these  comparisons. 
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where 

fr  is  the  bidirectional -reflectance  distribution  function 
(BRDF),10 

P0  is  the  incident  power. 

P  is  the  scattered  power, 
ft  is  the  solid  angle. 
k  *  2ir/X. 

X  is  the  wavelength  of  the  radiation, 

F  is  the  optical  factor  (given  explicitly  below), 

W{p)  is  the  power  or  the  variance  spectrum  of  the  surface 
height  distribution, 

p  *  It  (sin  0,  -  sin  0,)  is  the  surface  spatial  frequency, 
0,  is  the  angle  of  scatter, 

0,  is  the  incidence  or  specular  reflection  angle  with  re¬ 
spect  to  the  average  surface  normal, 
ss  means  that  the  source  and  the  receiver  both  are 
polarized  with  the  electric  vector  normal  to  the 
plane  of  incidence,  and 
pp  is  the  equivalent  parallel  situation. 

The  vector  theory  predicts  that 


THEORETICAL  BACKGROUND 

Classical  scalar  theory  has  been  used  quite  successfully  to 
calculate  diffraction  from  obstacles  and  apertures.  As¬ 
sumptions  are  made  so  that  the  properties  of  the  materials 
(such  as  conductivity  and  dielectric  constants)  do  not  enter 
the  calculations.  No  notice  is  taken  of  polarization,  but  if 
there  exists  only  a  single  polarization  in  the  incident  beam  and 
ellipticity  is  not  introduced,  the  problem  is  a  scalar  one.  For 
the  metallic  sample  and  the  angles  of  incidence  that  are  used, 
vector  theory  shows  that  the  ellipticity  introduced  on  reflec¬ 
tion  is  negligible. 

Vector  theory,  using  only  the  first  term  of  the  scattered -field 
expression  as  described  by  the  references  cited,®  car  be  written 


Ftt  *  cos  0,  cos  0, ,  (3) 

-  sin  0,  sin  0,V* 

—  - - —  •  (4) 

cos  a,  cos  0,  / 

Scalar  theory  predicts,  for  instance,  from  Davies.1  that 

F  *  (cos  0,  +  cos  0,  )4/co$  0,.  (5) 


F„  (— 

I 


The  surface -height  spectrum  W(p)  is  the  Fourier  transform 
of  the  autocorrelation  of  the  surface  height  distribution.  This 
is,  of  course,  equivalent  to  the  spectrum  of  the  distribution  of 
the  variance  from  the  mean  surface  height.  Many  authors 
have  assumed  either  Gaussian  or  exponential  forms  for  the 
autocorrelation  function,  and  these  lead  to  Gaussian  or  Lo- 
rent2ian  power  spectra.  We  make  no  such  assumption  mi- 
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dally.  Simple  reasoning  leads  to  the  conclusion  that  a  two- 
dimensional  power  spectrum  should  be  used,  but  there  are 
reasons  for  concluding  that  a  one-dimensional  spectrum  is  also 
valid.  Appendix  A  shows  that,  in  the  latter  case,  A4  in  the 
BRDF  equations  should  be  replaced  by  A3. 

Since  we  measured  the  scattering  of  unpolarired  radiation 
with  an  unpolarized  receiver  as  well  as  the  ss  and  the  pp  sit¬ 
uations,  the  theoretical  relationships  that  we  used  should  be 
given  It  is  well  known  that  the  total  flux  is  given  by 

d Puu  *  -  (d P"  +  dP,p  +  d P*'  +  d Pop).  (6) 

2 

However,  when  polarized  radiation  is  obtained  from  un¬ 
polarized  radiation  with  a  polarizer  pair  that  has  a  transmis¬ 
sion  of  either  r„  or  rpp,  one  has 

d P"  *  r„dP“",  (7) 

dPpp  «  TppdPw.  (8) 

When  these  substitutions  are  made  in  Eq.  (6)  the  result  is 

-  l-  <J,“  +  /,"’>•  (9) 

We  have  also  made  use  of  the  fact  that  the  cross  terms  are 
theoretically  zero  and  we  found  them  to  be  experimentally 
unmeasurable. 

The  same  result  is  obtained  from  scalar  theory  but  from 
different  physical  reasoning.  Since  scalar  theory  treats  all 
polarizations  equally,  the  BRDF  is  always  the  average  of  all 
the  polarization  terms.  For  in-plane  measurements  the  result 
is  given  above.  However,  for  out-of-plane  measurements  the 
scalar  result  should  be  half  of  the  vector  result 

EXPERIMENTAL  PROCEDURE  AND  RESULTS 

All  data  were  taken  on  a  single  sample  of  well-polished  fused 
silica  that  had  a  70-nm -thick  evaporated  coating  of  aluminum 
on  its  surface.  Such  a  surface  has  often  been  called  super- 
polished.  The  sample  was  measured  on  the  Naval  Weapons 
Center  Talystep  instrument11  and  found  to  have  an  rms 
roughness  of  about  8  A  (0.8  nm.  or  8  X  10~4  mid).  The  scat- 
terometer  used  for  the  measurements  has  been  described 
elsewhere.12  Measurements  were  made  at  0.6328,  3.39.  and 
10  6  wm  at  angles  of  incidence  of  10®,  30°,  and  50°  for  ss  and 
su  polarization  and  0.488,  T5145,  and  0.6328  pm  for  pp  po- 


Fig  1  BRDF  d«u  for  0  6328 3.39  Mm(0).  and  10.6 

for  10*  incidence  and  the  corresponding  one-dimensional  vector 

theorv  u.o  <  solid  lines!. 
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Fig.  2.  BRDF  data  for  0.6328  urn  (*).  3.39  am  (0).  and  10.6  Mm  (  +  ) 
for  3Cr  incidence  and  corresponding  one-dimensional  vector  theory 
(ss)  < -olid  lines). 


Fig.  3.  BRDF  data  for  0  6328  <*).  3.39  mo  (0).  and  10.6  Mm  (4-) 

for  50°  incidence  and  corresponding  one-dimensional  vector  theory 
1 55)  (solid  lines). 


Fig.  4.  BRDF  data  for  0.6328  Mm  ( +).  0.5145  Mm  (x).  and  0.4880  Mm 
(0).  and  the  corresponding  one-dimension&J  vector  theory  (pp)  (solid 
lines)  for  30*  incidence. 

larization.  The  perpendicular  polarizations  were  not  available 
at  the  longer  wavelengths  because  of  equipment  limitations. 
It  is  generally  true  that  the  data  are  reliable  from  0  -  0O 
(sin  6,  —  sin  0, )  values  greater  than  about  0.08  and  until  the 
data  curl  up  at  large  angles  because  of  noise  or  other  extra¬ 
neous  effects.  Figures  1-5  show  representative  data  and 
vector -theory  curves  that  are  adjusted  by  the  choice  of  rms 
roughness  and  autocorrelation  length;  Figs.  8-8  show  the  in¬ 
strument  function  and  the  noise  level.  We  believe  that  the 
uncertainty  of  our  measurements  is  20^t  in  most  cases,  al- 
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Fig.  5.  BRDF  for  0.632S-mU)  un polarized  theory  exponential  model 
for  50*  incidence 


Fig.  6-  Instrument  function  for  positive  I  +)  and  negative  (*)0  - 
at  0.6328  M*n  and  noise  level  for  positive  (0)  aiid  negative  (1)  d  - 

do ■ 


Fig.  7.  Instrument  function  for  positive  (  +  >  and  negative  (*)  3  -  do 
and  noise  level  for  positive  (0)  and  negative  (1)  3  -  do  at  3.39  pm. 

though  it  can  increase  when  we  approach  the  noise  at  low 
values  of  BRDF  and  at  high  or  low  values  of  J  -  do- 

RESULTS 

One  of  the  most  satisfying  ways  of  comparing  the  data  and  the 
theory  is  to  calculate  the  surface  height  spectrum  from  the 
BRDF  data  and  the  theoretical  relationship.  In  general,  this 
means  that 

dot 
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tin 

The  radian  wave  number  is  raised  to  either  the  third  or  the 
fourth  power,  depending  on  the  assumption  of  a  one-dimen¬ 
sional  or  a  two-dimensional  power  spectrum.  The  factor  F 
is  taken  from  the  appropriate  theory.  Nothing  more  is  as¬ 
sumed  about  either  W(p)  or  the  data. 


Fig.  8.  Instrument  function  for  positive  (+>  and  negative  (*)  d  -  do 
and  noise  ievei  for  negative  3  —  do  at  10.6  urn. 


Fig.  9.  Power  spectrum  for  0.6328. 3.39,  and  10.6  urn  at  30*  incidence 
with  two-dimensional  Beckmann  model. 


Fig.  10.  Power  spectrum  for  0.6328.  3.39.  and  10.6  nm  at  30*  inci¬ 
dence  with  one -dimensional  Beckmann  model. 
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Fig.  11  Power  spectrum  for  0.6328.  3.39.  and  10.6  mo  at  30*  inci¬ 
dence  with  two-dimensional  Davies  model. 


Fig.  12.  Power  spectrum  for  0.6328.  3.39.  and  10.6  um  at  30*  inci¬ 
dence  with  one-dimensional  Davies  model 


Figure  9  shows  the  result  of  calculating  W2(p)  using  a 
two-dimensional  Beckmann  model.  Three  different  spectra 
for  the  actual  (unique)  surface- height  distribution  are  ob¬ 
tained.  This  is  clearly  unsatisfactory.  Figure  10  shows  a 
similar  situation  for  the  Beckmann  one-dimensional  model. 
The  Davies  model  preceded  that  of  Beckmann.  Both  are 
scalar  and  both  are  unsatisfactory,  as  Figs.  11  and  12  dem¬ 
onstrate. 

Figure  13  shows  the  situation  for  the  vector  theory  with  a 
two-dimensional  power  spectrum.  Although  it  is  a  better  fit, 
the  three  curves  are  different  by  about  an  order  of  magnitude 
in  the  central  portion  of  the  curve.  A  typical  experimental 
uncertainty  is  about  the  size  of  the  markers  (0.  * .  +■ )  showing 
the  data  points.  Figure  14  is  very  satisfying.  The  lO-jim  data 
are  indicated  by  a  +  and  extend  almost  to  p  *  1  jim In  that 
region  they  cur)  up.  This  is  due  to  experimental  noise  effects. 
Similar  but  less-pronounced  effects  occur  for  3.39-mh)  data 
a  little  beyond  p  *  1  and  for  0.6328-am  data  at  about  p  *  7. 
In  the  region  from  p  *  0.1  vm~]  top  *  1  urn"1,  where  all  the 
curves  overlap,  they  all  predict  the  same  surface-height  dis¬ 
tribution — a  straight  line  with  a  slope  that  is  approximately 
-2.  The  0.6328-Mm  data  extend  the  same  line  about  another 
order  of  magnitude  in  p.  The  other  end  of  the  curve  has  some 
difficulties  associated  with  it.  The  0.6328-  and  3.39-Mm  data 
seem  to  extend  to  shorter  values  of  p  with  the  same  slope. 


whereas  the  10-jjm  data  tend  to  roll  over.  Unfortunately  this 
is  the  region  in  which  the  measured  data  are  contaminated 
with  scattering  from  the  instrument  near  the  specular  direc¬ 
tion.  Note  that  p  ■  0.5  nm~l  corresponds  to  0  -  0q  *  0.05, 
0.27,  and  0.84  for  0.6328,  3.39,  and  10.6  pm.  respectively. 
Instrument  scatter  generally  starts  at  about  0  -  0O  «  0.07  (and 
influences  measurements  at  shorter  values),  but  the  influence 
on  the  measurement  is  a  function  of  the  absolute  value  of  the 
measurement.  A  theoretical  curve  is  also  shown  in  Fig.  14. 
It  is  a  one-dimensional  Lorentzian  derived  from  an  expo¬ 
nential  correlation  function  and  fitted  to  our  data.  The  rms 
surface  height  a  and  the  correlation  length  /  were  found  to  be 
33.5  A  and  10  nm,  respectively.  This  theoretical  curve  fits  our 
reliable  data  well,  but  other  spectra  could  be  used. 

Figures  15-19  are  similar  curves  for  other  angles  of  inci¬ 
dence.  Some  of  the  fits  are  not  so  good,  but  we  believe  that 
they  all  support  the  one-dimensional  vector  theory  quite  well. 
Figure  20  shows  similar  data  for  pp  polarization  but  for  three 
shorter  (and  more  similar)  wavelengths,  and  Fig.  21  shows  that 
of  uu  i  un)  polarization.  Both  are  very  close  to  the 
0.6328-um  power  spectrum.  The  data  support  the  conclu¬ 
sions  earlier  stated. 

Another  approach  was  taken  to  assess  the  validity  of  the 
Lorentzian  spectrum.  We  assumed  that  Wx (p )  was  Lorent- 


Fig.  13.  Power  spectrum  for  0.6328.  3  39,  and  10.6  um  for  30*  inci¬ 
dence  with  two-dimensional  vector  theory  and  corresponding  theo¬ 
retical  curve  (solid  line). 


Fig.  1 4.  Power  spectrum  for  0.6328. 3.39.  and  10.6  am  for  -30*  in¬ 
cidence  with  one-dimensional  vector  theory  and  corresponding  the¬ 
oretical  curve  (solid  line). 
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Fig.  15.  Power  spectrum  for  0.6328.  3  39.  and  10.6  Mm  for  - 10°  in¬ 
cidence  with  one-dimensional  vector  theory  and  corresponding  the¬ 
oretical  curve  (solid  line). 


Fig.  18.  Power  spectrum  for  0.6328.  3.39.  and  10.6  Mm  for  30*  inci¬ 
dence  with  one-dimensional  vector  theory  and  corresponding  theo¬ 
retical  curve  (solid  line). 


Fig.  16.  Power  spectrum  for  0.6328.  3.39.  and  10.6  un  for  -50*  in¬ 
cidence  with  one-dimensional  vector  theory  and  corresponding  the¬ 
oretical  curve  (solid  line). 


Fig.  19.  Power  spectrum  for  0.6328.  3.39.  and  10.6  Mm  for  50*  inci¬ 
dence  with  one-dimensional  vector  theory  and  corresponding  theo¬ 
retical  curve  (solid  line). 


Fig.  17.  Power  spectrum  for  0  6328. 3.39,  and  10.6  Mm  for  10*  inci¬ 
dence  with  une-dimensional  vector  theory  and  corresponding  theo¬ 
retical  curve  (solid  liner 

zian.  and  we  calculated  <i  and  /.  Table  l  shows  these  results 
for  n  obtained  by  using  vector  theory.  The  correlation  length 
was  found  to  be  constant  at  1 2  jim.  For  ss  data  three  values 


Fig.  20.  Power  spectrum  for  0.488.  0.5145.  and  0.6328  Mm  at  30* 
incidence  with  one-dimensional  vector  model  and  corresponding 
theoretical  curves. 

for  <j  were  obtained:  27.  40.  and  50  A.  For  pp  data  three 
values  were  found:  13. 19.  and  25  A.  For  uu  data  three  values 
were  found:  11.  15.  and  16  A.  These  were  all  found  by  curve 
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Fig-  21.  Power  spectrum  for  uu  polarization  at  30*  incidence  baaed 
on  one-dimensional  vector  model. 

fitting.  Table  2  shows  similar  results  when  scalar  theory  is 
used. 

The  pp  values  are  about  half  of  those  obtained  from  ss  data 
and  vector  theory.  We  are  told  that  this  is  not  unusual.  The 
uu  data  agree  well  with  the  pp  and  with  the  scalar- theory 
result*. 

Table  3  shows  results  obtained  by  curve  fitting  the  derived 
W(j3)  results  and  by  integrating  under  these  curves.  We 
would  expect  the  a  values  for  a  given  angle  of  incidence  to  be 
the  same  whether  they  were  obtained  by  Fitting  BRDF  or 
surface  spectrum  data.  The  integrals  are  obtained  by  adding 
the  area  under  the  power  spectrum  piece  by  piece.  Some 
differences  should  arise  from  this  approximation. 

The  rau  surface  height  obtained  for  ss  vector  theory  is 
about  3-1  nm  When  obtained  by  using  pp  data  it  is  about 
2  nm.  When  obtained  by  using  scalar  theory  it  is  about  1.5 
nm.  Independent  measurements  by  Bennett 1 1  using  a  surface 
profilometer  give  0.8  nm.  as  do  the  results  of  Wyant  using  a 
heterodyne  interference  microscope.'3 
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DISCUSSION 

It  should  be  emphasized  that  all  these  result*  are  for  an  ex¬ 
haustively  measured  single  sample. 

The  series  of  curves  of  W(p)  versus  p  shows  clearly  that 
only  the  one-dimensional  vector  theory  it  a  reasonable  model 
for  scattering  regardless  of  the  form  of  W(p).  There  is  at 
present  no  explanation  for  the  discrepancy  between  pp  and 
: .!  data  in  predicting  o,  although  our  pp  data  were  limited  in 
wavelength  range.  The  agreement,  or  lack  of  it.  among  the 
several  different  methods  of  obtaining  a  is  disturbing.  It 
seems  to  indicate  that  the  connection  between  the  BRDF  and 
the  surface  properties  ha*  not  been  completely  established. 

The  set  of  W^p)  curves  using  one-dimensional  vector 
theory  provides  encouraging  evidence  for  wavelength  and 
angle  scaling  for  values  of  p  larger  than  1//.  in  spite  of  the  in¬ 
complete  theoretical  verification. 

The  next  steps  for  investigation  are  pp  data  for  a  largtr- 
wavelength  region  and  more-reliable  measurements  for 
smaller  values  of  3  -  do.  It  seems  obvious  also  that  other 
samples  should  be  investigated. 

APPENDIX  A.  TRANSITION  FROM  TWO- 
DIMENSIONAL  TO  ONE-DIMENSIONAL 
POWER-SPECTRA  EXPRESSIONS 

From  Eq.  (1),  the  two-dimensional  BRDF  can  be  expressed 
as 

/,*■  *  BRDF  W2(p)  IA1) 

ir2 

The  power  spectrum  W2(p)  is  the  Fourier  transform  of  au¬ 
tocorrelation  function  g(r).  If  the  surface  is  assumed  to  be 
isotropic,  and  the  exponential  autocorrelation  function  is  also 
assumed,  the  explicit  form  of  Eq.  ill  is 

fr"  ’~ik‘x  F“  f.  f  .  exp(— fx.  > )//) 

x  exp(—;2x(x/i  +y/v))drdv.  (A2) 


Table  1.  Values  of  rau  Surface  Height  for  Different  Measurements  Obtained  with  Vector  Theory  (in  angstroms) 


Incidence 

Angie  (deg) 

Wavelength 

S3 

PP 

uu 

0.6328 

0.6328 

3.39 

10.6 

0.6328 

0.5145 

0.4880 

“10 

27 

41 

41 

25 

25 

25 

15 

*10 

°7 

27 

25 

25 

25 

“30 

41 

41 

41 

19 

19 

19 

16 

30 

41 

41 

41 

19 

19 

19 

11 

-50 

27 

50 

41 

13 

13 

13 

50 

o  - 

27 

27 

19 

19 

19 

Table  2.  Values  of  rms  Surface  Height  Obtained  by 
Fitting  Scalar  Theory  to  BRDF  Data 


Incidence 

Angle 

'deg) 

Wavelength 

ss 

PP 

uu 

0.6328 

06328 

3.39 

10.6 

0.6328 

3.39 

106 

“10 

8 

12 

10 

8 

9 

9 

8 

10 

9 

8 

8 

8 

9 

9 

5 

-30 

11 

15 

13 

8 

10 

9 

7 

30 

13 

13 

11 

10 

11 

11 

8 

-50 

15 

16 

13 

13 

15 

16 

11 

50 

16 

18 

21 

18 

16 

18 

11 

Table  3.  Values  of  rms  Surface  Heights  Obtained 
from  the  Vector  Theory  by  (A)  Integration  and  <B) 
Curve  Fitting  to  W{p) 


incidence 

Angle « deg) 

A 

B 

-10 

21 

28 

10 

IT 

23 

-30 

24 

34 

30 

20 

31 

-50 

32 

35 

50 

32 

43 

51 
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where  exp(-(x,  y)H\  is  the  two-dimensional  exponential  au¬ 
tocorrelation  function.  If  it  is  separable,  exp(-(x,  y)ll\  is 
equal  to  (expl-x/1)  exp(y/l)|. 

Next,  we  assume  that  the  surface  is  isotropic  in  one  di¬ 
mension  only,  the  j  direction.  We  can  integrate  over  d«  to 
obtain  the  form  of  dp/dd.  We  know  that  dll  is 


dQ  «  dd.do,  sin  d,  ■ 


(A3) 


We  have 


2  t  2t 

d q  *  2xd/y  *  —  dky  *  —  s*n  008  '  (A4) 


When  Eq.  (A4)  is  substituted  into  Eq.  (A3),  one  has 
dfl  sin  fl,d q  a  dq 


dO  s 


«  dd  x  ■ 


2ir 

—  sin  cos 


2 x 

—  cos  0t 


(AS) 


With  the  assumptions  of  separability  and  one-dimensional 
isotropicity.  dq/d0  can  be  obtained  by  integrating  Eq.  (A2) 
over  all  ranges  of  q  instead  of  <p  by  using  Eq.  (A5). 

Thus  Eq.  (A2)  becomes 


dd 


xxx: 


expl-x/I) 


—  cos  <o. 


But 


x  expl-;2irx/,)dx  expi -y/Oexpl -;2Ti/y )dy.  (A6) 

21 


x: 


exp(  —  e/1  )exp(  — ;  2irx/K  )dx  : 


x; 


exp(>7)exp(-./2;r\X>  )dy 


1  (2t/,I)j 

2 1 

1 1  +  (pt)2 ' 

21 

'  i  +  (2  x/yn3 
21 

1  +  1 qt )2‘ 


IA7) 


(A8) 


Equation  (A6)  becomes 


cLP  I6ir-  <r-  21  C~  21 

d«  *  a 4  0  i  *  (pir-  J--i  +  {qi)-dq' 


—  COS  d. 


However. 


f  *  - - - ;dq  -  f'-SL; 

J  -•  l  -t-  iq/i-  J--1  +  X- 


(let  ql  »  r). 

IA10> 


Equation  i  A9l  reduces  to 

d P  16*’ ,  „  ,  <?• 


'  -b 

A  1  COS  <$, 

.IgEfiri  2*lL- 

^‘Jl  1  4-  |p/)* 


2/ 


l  +  ip()J 

2  a:l 


1  +  (p/V- 

(cos  <p,  *  1  in  the  plane). 

(All) 
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Abstract 

Vector  theory  has  been  used  to  obtain  the  power  spectral  density  of  the  surface  height 
distribution  at  several  different  wavelengths,  conditions  of  polarization  and  geometries. 
The  surface  height  was  then  determined  both  by  integration  and  by  curve  fitting. 

Introduction 


One  of  the  goals  of  the  theory  of  surface  scattering  is  to  relate  the  scattered  light 
to  the  surface  properties.  Determining  both  the  correlation  length  2  and  the  rms  roughness 
3  of  a  random  surface  might  be  considered  the  main  problem  of  inverse  scattering,  although 
determining  scaling  laws  and  surface  distribution  functions  are  certainly  of  equal  concern. 

In  previous  publications  we  have  shown  rather  good  agreement  between  measurements  at 
several  wavelengths  for  different  incidence  angles  and  for  different  polarization  and  the 
vector  version  of  electromagnetic  theory. 1  Here  we  apply  it  to  the  determination  of  o  and 
f.  The  BRDF  written  as  f  can  be  written2 

fr  *  F:.  Fo  Fs  (1) 

F^  »  wavelength  factor 

Fq  m  optical  factor 

Fs  ■  surface  factor 

Appropriate  expressions  exist  for  each  one  of  these  factors.  The  optical  factor  is  a 
complicated  function  of  the  angles  of  incidence  and  scattering,  the  form  of  which  depends 
upon  the  polarization  conditions.  The  wavelength  factor  is  k3  or  k1*  depending  upon  whether 
the  surface  factor  is  assumed  to  be  a  one-  or  a  two-dimensional  power  spectrum  of  the  sur¬ 
face  height  distribution.  The  surface  height  power  spectrum  is  then  given  by  either  of  the 
two  expressions 


w:  (p) 

2  J2i 

l  +•  p2  2  2 

w; ; p  > 

2- j2  22 

(2) 

( i-p2  2  2  )  2/2 

The  subscript  indicates  whether  a  one-  or  a  two-dimensional  form  is  assumed. 

The  value  of  j  can  then  be  determined  by  several  methods:  a  fit  to  the  BRDF  data,  a  fit 
to  the  inferred  power  spectrum,  integration  under  the  asymptotic  approximation  or  a  two- 
waveie.ngth  ratio  technique.  We  discuss  each  technique  in  turn  and  then  compare  them. 

BRDF  curve  fitting 

We  have  taken  our  measured  data  in  the  past,  assumed  one  of  the  two  forms  (actually  W, 
worked  better)  of  the  surface  power  spectrum  and  used  trial  and  error  for  determining  a 
and  i.  This  result  should  be  expected  from  the  form  of  the  equations.  We  note  that  the 
determination  depends  upon  the  model  for  W. 

Power  spectrum  curve  fitting 

The  BRDF  data  can  be  divided  by  the  appropriate  surface  factor  to  obtain  the  power  spec¬ 
trum  of  the  surface.  We  did  this  for  several  different  wavelengths,  polarizations  and 
incidence  angles.  The  results  showed  that,  at  least  for  the  sample  we  used,  the  power  spec¬ 
trum  was  one-dimensional.  One  can  then  find  the  values  of  i  and  2  that  give  these  power 
spectrum  curves.  Alternately  one  can  integrate  under  the  curve  since  the  total  integral 
of  the  power  spectrum  is  .  Numerical  integration  is  tedious,  so  integration  under  the 


asymptotic  approximation,  described  next,  is  appropriate.  Note  that  with  numerical  inte¬ 
gration  or  curve  fitting  to  the  power  spectrum  only  the  dimensionality  of  W  is  assumed 
and  not  its  form. 


Asymptotic  methods 


If  the  power  spectrum  can  be  properly  represented  by  the  one-dimensional  model,  as  our 
sample  can  be,  the  curve  can  be  represented  asymptotically  by  two  straight  lines  that 
intersect  at  (p;,h).  The  integral,  as  shown  in  Appendix  A  is  2hpi  and  this  is  equal  to 
2 tu2.  We  note  that  if  the  spectrum  is  two-dimensional  the  slope  is  -3  rather  than  -2  and 
the  result  is  different  but  easily  determined.  This  method  depends  upon  both  the  dimension¬ 
ality  and  the  model. 

Ratio  technique  for  calculating  1 


We  now  consider  determining  a  based  on  two  measurements  at  the  same  value  of  3-3  but 
at  different  wavelengths.  It  can  be  seen  immediately  that  the  surface  factors  drop’out. 
For  either  of  the  models  we  discussed  3  also  drops  out.  The  expressions  have  only  X  and 
3-3  in  them.  Either  complicated  algebraic  manipulations  or  numerical  techniques  can  then 
be  used  to  obtain  l. 


Table  1  gives  values  of  a  obtained  by  the  BRDF  fitting  technique  for  a  variety  of  angles, 
polarizations  and  wavelengths.  The  values  obtained  range  from  11  to  42,  a  factor  of 
almost  4.  This  is  probably  due  to  variations  in  individual  measurements  and  inconsistency 
in  the  (tedious)  curve  fitting.  Table  2  gives  a  values  obtained  by  three  other  different 
methods  for  the  best  estimate  of  the  power  spectrum  using  all  three  wavelengths.  The 
agreement  is  better,  but  we  should  recognize  the  smoothing  that  results  from  using  the 
three  wavelengths.  We  also  should  note  that  no  pp  or  uu  data  are  included.  Table  3  gives 
the  correlation  length  obtained  by  the  ratio  technique.  Table  4  gives  values  of  the  cor¬ 
relation  length  obtained  by  four  different  methods.  Table  5  gives  the  rms  surface  height 
obtained  by  eight  different  ways.  The  last  two  are  independent  measurements  made  by  Jim 
Wyant  with  a  heterodyne  interferometer  and  by  Jean  Bennett  using  a  mechanical  probe. 


Table  1.  RMS  surface  height  found  from  curve  fittin 


Angle  of 
incidence 


Table  2.  RMS  surface  height  (A) 


Area  integration  Curve-fitting  on 

of  oower  soectrum  (ss)  ss  oower  soectrum 


Asymptotic  area  integration 
of  oower  soectrum 


54 


Table  3.  Correlation  length  from  ratio  technique 

Mum) 


0.6328  —  3.39 

10* 

9.09 

3.39  - -  10.6 

10° 

10.08 

0.6328  < — -  10.6 

10« 

9.94 

3.39  -  10.6 

-10° 

8.55 

0.6328  - — •  10.6 

-10° 

10.74 

0.6328  - — -  3.39 

-30" 

8.91 

Average 

9.55 

Table 

4. 

Correlation  length  by 

four  methods 

Average  of  i 
by  ratio  technique 

9.55 

Asymptotic  method 
at  long  wavelength 

12 

Curve-fitting 
on  as  BRDF 

10 

Curve-fitting 

on  ss  power  spectrum 

10 

Table  5. 

Determination  of  rms  surf. 

ace  height 

o(A) 

area  under  ss  power  spectrum 

24 . 4 

asymptotic  area  under  ss 
power  spectrum 

26.4 

curve-fitting  of 
ss  BRDF  data 

32.3 

curve-fitting  of 
ss  power  spectrum 

34 . 9 

curve  fitting  of  op 
power  spectrum 

19  ' 

curve- fitting  of  uu 
power  spectrum 

14.1 

Wyant 

3 

Bennett 

7.5 

Discussion 


The  autocorrelation  length  can  be  found  a  number  of  different  ways.  Some  are  dependent 
upon  the  entire  model,  some  are  independent  of  the  surface  factor  (as  long  as  it  is  inde¬ 
pendent  of  wavelength) .  All  methods  gave  good  agreement.  They  also  agreed  with  an  inde¬ 
pendent  optical  method,  but  they  disagreed  considerably  with  a  mechanical  method.  The 
disagreement  is  still  a  matter  of  further  study. 

Several  methods  are  available  for  obtaining  the  rms  height.  The  method  of  BRDF  curve 
fitting  i3  tedious  and  completely  model  dependent.  Power  spectrum  fitting  requires  no 
assumption  about  the  surface  factor.  An  assumption  that  the  autocorrelation  is  exponen¬ 
tially  decreasing  is  necessary  for  either  fitting  the  curves  or  their  asymptotes.  The 
latter  method  is  easy  and  seems  as  accurate  as  the  data  warrant.  The  ss* data  give  a  values 
consistently  higher  than  the  other  methods.  This  disagreement  has  been  noted  before  and 
is  still  a  matter  of  conjecture. 

Table  6  provides  a  summary  of  the  methods  for  determining  the  autocorrelation  length, 
the  source  of  data,  the  required  assumptions  and  the  apparent  limitations.  Table  7  Is  a 


similar  presentation  for 

the 

rms  surface  height 

Table 

6. 

Methods  for  determining  correlation  length 

• 

3RDF  curve  fitting 

BRDF  data 

Optical  factor 
surface  factor 

Curve  fitting  error 

Power  spectrum 
curve-fitting 

Power  spectrum 

Optical  factor 

Curve  fitting  error 

• 

Asymptotic  in 
long  wavelength 

Power  spectrum 

Optical  factor 

Instrumental 

scattering 

Ratio  technique 

BRDF  data 

Surface  factor 

Measurement  error. 

Limited  available 
data 

Table  7 

Methods  for  determining  surface 

height 

Data  source 

Assumption 

Limitations 

• 

BRDF  curve  fitting 

BRDF  data 

Optical  factor 

Curve-fitting  error 

Power  spectrum 
curve- fitting 

Power  spectrum 

Optical  factor 

Curve-fitting  erro- 

• 

Area  integration  of 
power  spectrum 

Power  spectrum 

Optical  factor 

Instrumental  Scatter¬ 
ing  time  consuming: 
tedious 

Asymptotic  area 
of  power  spectrum 

Power  spectrum 

Optical  factor 
Surface  factor 

Instrumental  Scatter- 
inn  :  fast,  easv 

Appendix  A.  Asymptotic  integration 

When  the  surface  spectrum  can  be  represented  by  two  straight  lines  on  a  log-log  plot, 
the  integration  is  relatively  simple.  The  area  under  the  curve  is  given  by 


where  K 
m 

rhis  reduces  to 


j  h  dp 
'o 


constant 
2  or  3 


-P: 

i  Z*  dp 

P:  P 


1 1  -  h  p •  *  K  -  —\ 

1  \P:  PzJ 

4  *  h  +  I  K  (j?'  j?) 
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for  nose  cases  (all  cases  we  have  considered)  p2  is  much  greater  than  pt  and  the  second 
term  of  the  parenthesis  may  be  ignored.  At  the  intersection  point 


h  * 

K/p i 2  or  h  » 

K/p : 3 

(6) 

hp  ■  m  K/p  !  hp  i  * 

Thus  the  two  areas  reduce  to 

K/p  ,  2 

(7) 

2  h  pi 

(8) 

A2  * 

1.5  h  p. 

(9) 

These  areas  are 

equal  to  2ro2 

so  that 

0  i  » 

(h  p!/'r)‘s 

(10) 

o2  * 

(3  h  p  x /4m )  ** 

(11) 

Appendix  3.  Ratio  technique  for  calculating  1 

We  assume  that  BRDF  measurements  have  been  made  at  each  of  two  wavelengths.  Then 

f  »  -  k !  -  F  2  a2i/[l+k;2  ( S-3  )  2  i 2  ]  (12) 

r  l  T  o  0 

f  =>  |  k23  Fq  2c2i/(l+k23(8-ao)Ji2)  (13) 

The  ratio  R  is  then  found  to  be 


R  * 

3 

\  u  +  *22 

! 3-3„ ) 2 i2 1 

(14) 

\»l  J 

'  tl  *  k i 2 

(8-3  )222] 

0 

Let  R, 

A 

i  2/x 

1 

Then 

R 

R,  3 

A 

1  t  k 

'.2(a-3o>2] 

«  ^1  +  k22(3-3o)2j 

(15) 

r  2  0  n2“| 

R 

R.  2 

1  » 

(2ml  2 (3-3 

0 

)2|(w  '  (w  1 12 

(16) 

(»A  - 1) 

(i I2- A  22 ) 

j  2 

(17) 

12t(3-3o)J 

2 (A i2-»22) 

Although  further 

manipulation  is 

possible,  the  solution  for  2  is  apparent. 

It  is  given  in 

terms  of  the  two  known  wavelengths  and  the  value  of  3-3  .  A  similar  result  is  obtainable 
for  the  two-dimensional  case. 
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